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We analyze solutions to loop-truncated Schwinger-Dyson equations in massless 

J\f = 2 and AA = 4 Wess-Zumino matrix quantum mechanics at finite temperature, 

<^ where conventional perturbation theory breaks down due to IR divergences. We find 

_J a rather intricate low temperature expansion that involves fractional power scaling in 

the temperature, based on a consistent "soft collinear" approximation. We conjecture 

that at least in the AA = 4 matrix quantum mechanics, such scaling behavior holds 

to all perturbative orders in the 1/N expansion. We discuss some preliminary results 

^ in analyzing the gauged supersymmetric quantum mechanics using Schwinger-Dyson 

equations, and comment on the connection to metastable microstates of black holes in 

the holographic dual of BFSS matrix quantum mechanics. 
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1 Introduction 

The sixteen supercharge SU{N) gauged matrix quantum mechanics was famously proposed 
by Banks, Fischler, Shenker and Susskind to give a non-perturbative formulation of M- 
theory [T] . It was later recognized [21 El IH E] that the precise holographic dual of the matrix 
quantum mechanics is an asymptotically null compactification of M-theory along the general 
lines of gauge/gravity duality [El [71 [8] . Despite intense efforts over a number of years (see 
f9] and references therein), most of the successful tests of the duality were realized to be 
consequences of supersymmetric non-renormalization theorems fiUi [TT| [T2| [T3l [H], while 
unsuccessful attempts [ISl [IS [IZj in testing the duality involved comparisons of the gravity 
and gauge theory results in regimes that do not overlap. It had become clear that in order 
to explore the semi-classical gravity (either type IIA or 1 1-dimensional supergravity) regime 
in the bulk, one must work in the genuinely strong coupling regime of the matrix quantum 
mechanics. 

A strong coupling test of the duality within the 't Hooft scaling regime was performed by 
numerically computing the free energy of the matrix quantum mechanics at finite tempera- 
ture using Monte Carlo method [IHl UHl EOl Ell [22] ■ The result was shown to be consistent 
with the expected free energy of black holes in the gravity dual. In order for the black hole 
horizon to lie in the semi-classical gravity regime, one needs to take the low temperature 
and large N limits of the matrix quantum mechanics. However, these are precisely the limits 
where the Monte Carlo computation becomes costly. 

An alternative approach to the thermal free energy of the matrix quantum mechanics 
was pioneered by Kabat, Lifschytz and Lowe, using the one-loop truncated Schwinger-Dyson 
equations (231 [211 [25] . This is the approach we will follow, and refine, in this paper. In place of 
the numerical solutions to the Schwinger-Dyson equations investigated in [221 1211 [25] , we will 
analyze the solutions to these equations analytically. We will see that in the low temperature 
limit, intricate cancelation occurs due to an approximate restoration of supersymmetry, 
which allows for nontrivial scaling behavior with temperature. Despite the fact that the 
naive 't Hooft coupling goes to infinity in the low temperature limit, the effective coupling 



parameter that controls the loop expansion of Schwinger-Dyson equation could be finite, and 
the scaling behavior of solutions to the one-loop Schwinger-Dyson equation could survive to 
all order in the 1/A^ expansion, when there are sufficiently many supersymmetries. 

Before describing our results, let us briefiy review the connection between black holes in 
the bulk geometry and the thermal free energy of the dual quantum mechanics. The gravity 
dual of the SU{N) BFSS matrix quantum mechanics at finite temperature is given by type 
IIA string theory in the near horizon limit of a near extremal black hole carrying N units of 
DO-brane charge [21 [3l H |5] 

ds^jj^ = -f-'/^Adf + f'/^A-\dr^ + r^dnl), 

_ CpgsNll 3 rl 
J = ^ — , Co = bOvr , A = l -, 

where dVL\ is the metric on a unit eight-sphere, and the black hole horizon is at r = tq. The 
Hawking temperature of the black hole is 

where Qym = 5's/4vr^/g. The free energy of the black hole is [26] 

/921q2p:7_14\ 1/5 / rp \ 9/5 

The black hole horizon lies in the type IIA supergravity regime when gy^N^^'^ <^ T <^ 
{dYM^Y^^- Working in the 't Hooft scaling limit on the gauged quantum mechanics side, 
one takes the large N limit while keeping the dimensionless 't Hooft coupling gyj^jN/T^ 
finite. The T^'^ scaling of the free energy or the entropy of the black hole is expected to hold 
in the matrix quantum mechanics in the limit of large 't Hooft coupling, or equivalently, the 
low temperature limit. 

One must be cautious about the meaning of the free energy of the matrix quantum 
mechanics at large A^. Since BFSS quantum mechanics has exactly fiat directions, at finite 
A^, the free energy is infinite, refiecting the continuum of scattering states. There are only 
A^ fiat directions, however, whereas the entropy of the black hole in the gravity dual scales 
like A^^ times a function of the dimensionless 't Hooft coupling. If we regularize the volume 
divergence by an IR cutoff at distances much greater than the horizon size of the black hole 
but much smaller than e^, we expect the black hole states to dominate the contribution to 
the free energy. In particular, if we take the infinite A^ limit first, and then take the volume 



to infinity while maintaining that the volume grows no faster than exponential in N, the A^^ 
coefficient of the free energy should remain a finite function of the 't Hooft coupling. This is 
the planar free energy of interest. The corresponding, finite, entropy of order N"^ generally 
counts metastable states rather than exact energy eigenstates in the quantum mechanics. 
Nonetheless, the lifetime of these metastable are expected to go to infinity (exponentially) 
in the infinity A^ limit. 

In the gravity dual, the metastable states are the microstates of the DO black hole, which 
decays by emitting DO-branes via Hawking radiation. Note that the DO-branes and their 
BPS bound states are the only particles that can escape to infinity in the type IIA 0-brane 
decoupling geometry. To see this, consider the Born-Infeld effective action for a probe DO- 



brane in the background (1.1) 

I dt f-^A^'^^l + fA-'^r'^- I dt f- 
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While this effective action a priori holds only in the supergravity regime, there is strong 
evidence that it in fact continues to hold at arbitrary large radial distance where the curvature 
is at string scale (while the string coupling goes to zero). This is presumably due to the 
supersymmetry preserved by the asymptotic geometry. Note that in the large r limit, the 
mass of the DO-brane cancels the potential energy, and the action reduces to that of a non- 
relativistic particle. The situation is in contrast to the Hawking radiation of black holes in 
global AdS where all emitted particles bounce back in finite time, making it possible for the 
black hole to be in equilibrium with a thermal bath. We conclude that the black hole in the 
0-brane decoupling geometry is unstable via the emission of DO-branes. The Hawking decay 
rate of emitting a DO-brane is computed in Appendix [Aj One finds 

( 359 2^ / 2r9 

IpJ 

where Ip = {2'KgsY^^ls is the IID Planck length and i?io = Qsh is the radius of the M-theory 
circle. The exponential factor can be understood in terms of the chemical potential for the 
DO-brane charge. This formula is valid in the type IIA supergravity as well as in the M- 
theory regime where the lift of the black hole solution to 11 dimensions is thermodynamically 
stable. The latter is expected when tq ^ N^/^lp. This is also the regime where the decay 
rate is exponentially suppressed. 

The metastable microstates of the DO black hole should be counted, to leading order 
in the 1/A^ expansion, by the planar free energy of the matrix quantum mechanics. Even 
in the high temperature regime, where the black hole horizon spills into the stringy part 
of the bulk geometry, and where one naively expects the matrix quantum mechanics to 



be weakly coupled, one encounters infrared divergences in conventional perturbation theory. 
Similar IR divergences were previously encountered in the two-loop computation of scattering 
amplitudes on the Coulomb branch of the theory, and have been essentially ignored ^7\ [2H] • 
The IR divergence can be cured (non-perturbatively) if one solves for the exact propagators 
using Schwinger-Dyson equations. This is the approach of [23l IMl ES], where the authors 
studied the self-energies of various fields, as well as the free energy of the theory at finite 
temperature, using the one-loop truncated Schwinger-Dyson equations, which amounts to 
a mean-field approximation. In |23l [211 125], the solutions of the self-energies were found 
numerically. While in a certain temperature range the result seemed to be consistent with 
the expectation from the gravity side, the numerical solution appears to break down below 



a certain temperature. From the gravity side, it is clear that (1.3) should be valid for 
T /{gyM^Y^^ ^ A^"^/^^, and thus the temperature can be taken to be arbitrarily small in 
the 't Hooft limit. It would be highly desirable to have an analytic understanding of the 
T^/^ scaling of the low temperature free energy /entropy in the gauged quantum mechanics. 

Though broken at finite temperature, supersymmetry plays an important role in the 
low temperature limit of the solutions to the Schwinger-Dyson equations. As pointed out in 
[23l EH |25] , in working with a truncated set of Schwinger-Dyson equations, either to a certain 
loop order or by including a certain finite subset of renormalized vertices, the equations must 
be manifestly supersymmetric in the zero temperature limit (though the solutions may be 
singular in the zero temperature limit) in order to have any chance of capturing the correct 
low temperature limiting behavior. In particular, the self-energies of the auxiliary fields in a 
supermultiplet must be included in the S-D equations. In a gauge theory, the solution to a 
truncated set of S-D equations will depend on the choice of gauge. While the solution to the 
exact S-D equations should clearly be independent of the gauge-fixing condition, in working 
with the truncated S-D equations, a manifestly (off-shell) supersymmetric gauge-fixing is 
necessary. The familiar Wess-Zumino gauge breaks all supersymmetries, and cannot be 
applied for our purpose, namely to extract the low temperature physics from the solutions to 
the truncated S-D equations. The authors of [231 1211 125] considered a gauge fixing condition 
that preserves manifest Af = 2 supersymmetries, out of the Af = 16 supersymmetries of 
BFSS matrix quantum mechanics. This leads to rather unconventional kinetic terms and 
fermion coupling. Alternatively, one may choose to work with gauge fixing conditions that 
manifest A/" = 4 or A/" = 8 supersymmetries. As a matter of fact, the A/" = 4 gauge fixing 
results in a rather complicated looking, fully nonlinear, action, and the A/" = 8 gauge fixing 
based on harmonic superspace requires the inclusion of infinitely many auxiliary fields. 

As a first step towards understanding the low temperature scaling behavior in BFSS 
matrix quantum mechanics, we study its truncation to the matter sector in the cases of A^ = 2 
and A/" = 4 gauge fixing. These may also be thought of as supersymmetric deformations of the 
BFSS matrix quantum mechanics, by turning off the coupling of the Af = 2 or Af = 4 gauge 

6 



multiplet. For instance, the truncation to A/" = 4 matter multiplets results in the A/" = 4 
quantum mechanics with three matrix matter multiplets $" and the cubic superpotential 
W = — -^e^^'^Tr ($"$^<l>'^). We refer to such a theory as an A/" = 4 Wess-Zumino (matrix) 
quantum mechanics. The Schwinger-Dyson equations of the A/" = 4 Wess-Zumino quantum 
mechanics have the particularly nice property that, going beyond the one- loop truncation, 
there are in fact no two-loop contributions, and there are no three-loop planar contributions 
(and the first planar correction to the one-loop S-D equations shows up at four-loop order). 
It is conceivable that the solution to the one-loop S-D equations in fact captures the correct 
scaling behavior of the planar free energy. 

In this paper, we will find an analytic low temperature expansion of the solution to 
the one-loop truncated Schwinger-Dyson equations for the A/" = 4 Wess-Zumino quantum 
mechanics. The key observation will be that loops containing the zero mode of the bosonic 
field (f)"' dominate the contribution to the self-energies of the nonzero modes and the auxiliary 
zero mode. This allows for the solution of the nonzero mode and auxiliary zero mode self- 
energies in terms of the boson zero mode self-energy. The S-D equation for the boson zero 
mode self-energy, on the other hand, is nontrivial only if we work to the next-to-next-to- 
leading order contributions in the low temperature expansion. In the end we find a nontrivial 
scaling behavior of the boson zero mode self-energy, which controls the scaling of the self- 
energies of all other modes. 

The matter multiplet of A/" = 4 Wess-Zumino quantum mechanics contains the following 
component fields, schematically: the boson 0, the fermion ip, and the auxiliary field /. At 
finite temperature, i.e. in the Euclidean theory where the Euclidean time is compactified 
with periodicity /3 = 1/T, let the self-energies for the momentum modes of 0, / and ip on 
the Euclidean time circle be an, rjn, and hr, respectively. Here n is an integer, whereas r is a 
half integer, reflecting the anti-periodic thermal boundary condition for the fermionic field. 
The self-energies for the nonzero modes are solved in terms of the boson zero mode ao, with 
the following results: 



CTr 
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27rn| / 2 
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where C„ = ^^ ,Q^sign(A;)sign(n — k)/k. Here we expressed the results in units where the 
dimensionful 't Hooft coupling k'^N is set to 1. For the zero mode self-energies (Tq and tjq, of 



the boson and auxiliary field /, we will find 
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ao = 2(^)2/5/3-7/5 + 0(^-2) 



Somewhat surprisingly, the low temperature expansion parameter is a fractional power of the 
temperature, namely 13~'^/^. The planar free energy can then be computed in the mean-field 
approximation |23]. The result is 

I3F = const - -{-fl'^N^Nfp-^''' + C»(/3-^/^). (1.8) 

where Nf is the number of chiral superfields. Curiously, this result differs from the expected 
scaling of BFSS matrix quantum mechanics (1.3) by one power of the expansion parameter 

/3-3/5. 

While we have not evaluated explicitly the four-loop correction to the planar Schwinger- 
Dyson equations and to the free energy, it appears that such higher-loop corrections could 
contribute at the same order as the one-loop contributions to the self-energies, in the low 
temperature limit, despite the fact that the naive dimensionless 't Hooft coupling is infinite in 
this limit. It is conceivable that the T^/^ scaling is exact for the large N J\f = A Wess-Zumino 
matrix quantum mechanics with cubic superpotential, and that as we continuously deform 
the BFSS matrix quantum mechanics to the Af = 4 Wess-Zumino model by turning off the 
gauge coupling, the low temperature scaling behavior of the planar free energy interpolates 
between T^/^ and T^/^ 

One might be puzzled by the following. Consider Model I, the N' = 4 Wess-Zumino matrix 
quantum mechanics with a single matter multiplet $ and superpotential W = Tr$'^, and 
Model II, the A/" = 4 Wess-Zumino matrix quantum mechanics with three matter multiplets 
X,Y,Z and superpotential W = Tt{XYZ). In the former case the spectrum is gapped, 
whereas in the latter case there is a continuum of scattering states due to flat directions. 
At finite A^, the free energy of the former should be exponentially suppressed in the low 
temperature limit, whereas that of the latter diverges due to the continuous spectrum. How 
are these consistent with our claimed scaling as follows from the Schwinger-Dyson equations? 
Our result suggests that in Model I, even though the spectrum is gapped at finite A^, in the 
large A^ limit the gap in the spectrum becomes very small, and if we take A^ to infinity first 
and then take the low temperature (or equivalently, strong 't Hooft coupling) limit, the free 
energy exhibits power scaling in the temperature. In Model II, on the other hand, there 
are different solutions to the Schwinger-Dyson equation. The scaling solution we described 
above, in particular, treats X, Y, Z on equal footing. There are other, singular, solutions 
that sets the self-energy of one of X,Y ot Z to zero and giving infinite self-energies to the 



remaining two fields. These singular solutions describe the phase of the theory where one of 
the three fields acquires a large expectation value, while the other two fields are very massive, 
in contrast to the "unbroken phase" described by the scaling solution. We conjecture that 
these distinct phases exist in the infinite A^ limit, and the tunneling between different phases 
are exponentially suppressed in the large A^ limit. 

When vector multiplets, along with ghosts in the supersymmetric gauge fixing, are in- 
cluded in the one-loop truncated Schwinger-Dyson equations, it appears that our low tem- 
perature expansion scheme is spoiled. It is not clear to us whether the one-loop S-D equation 
captures the correct low temperature physics in this case. It is likely that loop corrections 
to the quartic and possibly higher vertices of the vector multiplet, which are not taken into 
account by the one-loop S-D equation, are needed to obtain a nontrivial low temperature 
scaling. We will nonetheless discuss preliminary results on the low temperature effective ac- 
tion of the vector multiplet by integrating out matter multiplets using the one-loop truncated 
S-D equations, as well as the high temperature expansion. The hope is that an improved 
S-D equation for vector multiplets will produce the T^/^ scaling of BFSS matrix quantum 
mechanics. Perhaps more auxiliary fields need to be included in the S-D equations, or a 
gauge fixing that preserves more manifest supersymmetries is needed. This is left for the 
future. 

The rest of the paper is organized as follows. In Section |2| we introduce the M = 2 
Wess-Zumino matrix quantum mechanics and derive its one-loop truncated Schwinger-Dyson 
equations. In Section|3} we solve the Schwinger-Dyson equations in the low temperature limit 
using a "soft collinear" approximation scheme, and compute the mean-field free energy. We 
then investigate the corrections from higher-loop diagrams, and argue that although the 
temperature scalings appear to be spoiled by these corrections in the A/" = 2 Wess-Zumino 
matrix quantum mechanics, they should remain valid in the M = A version of the theory 
in the planar limit. In Section |4| we introduce the A/" = 4 Wess-Zumino matrix quantum 
mechanics, and repeat the low temperature analysis. In Section |5| we discuss how our results 
for the Wess-Zumino quantum mechanics can be extended to the full BFSS matrix theory by 
coupling it to a vector multiplet; in particular, we present the supersymmetric gauge-fixing 
conditions and write down the corresponding Schwinger-Dyson equations. We also discuss 
the various phases of solutions. In Section [6} we explore the high temperature limit of BFSS. 
In Section [7| we discuss future prospects of this program, including ways to write down the 
BFSS action that preserves more manifest supersymmetries, and applications of our methods 
to supersymmetric quantum field theories in other dimensions such as the two-dimensional 
M = (2, 2) Landau- Ginzburg model. The Hawking decay rate of the black hole in the 0-brane 
decoupling geometry is derived in Appendix |Xj Details on the convention of A/" = 2 and A/" = 
4 superspace, the one-loop Schwinger-Dyson equations, and the low temperature expansion 
of the solutions are given in Appendices [BJJF} Finally, the high temperature expansion of the 
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Schwinger-Dyson equations for BFSS matrix theory is analyzed in Appendix [G 



2 M = 2 Wess-Zumino Matrix Quantum Mechanics 

In the ID A/" = 2 language, the BFSS matrix quantum mechanics consists of one vector 
multiplet and seven matter multiplets. Its truncation to the matter sector is an A^ = 2 
Wess-Zumino matrix quantum mechanics with the 7 matter multiplets interacting through 
a G'2-invariant cubic superpotential. The one-loop truncated Schwinger-Dyson equations for 
this system have been studied numerically in |23] . In this section we recall the form of the 
Schwinger-Dyson equations for the A/" = 2 Wess-Zumino model at finite temperature, and 
set up the notations for the analytic results in subsequent sections. 

2.1 The Action 

Consider M = 2 Wess-Zumino matrix quantum mechanics with flavor symmetry G. Let the 
matter multiplets be in some iVj- dimensional representation of G, labeled by index a. We 
will assume that G has a rank-3 totally antisymmetric invariant tensor e"'^'^ normalized by 



ubcubd „s:cd 



e 



^abd ^ ^^cd ^2.1) 



in this A^j-dimensional representation. If G is the 7 of G2 as for the matter sector of BFSS, 
then Nf = 7 and c = 3/2. The ID A/" = 2 superspace is introduced in Appendix Ib] In this 
language, the A/" = 2 Wess-Zumino quantum mechanics in Euclidean signature contains Nf 
real superfieldq^ 

^- = (f)- + z^-e^ + ire\ (2.2) 

interacting through a cubic superpotential 

W = -^e^^'^Tr (<l>n$^ $"]) . (2.3) 

6a/c 

Here $'' are matrix superfields in the adjoint representation of S't/(iV)n and in flavor sym- 
metry G. The coefficient for the superpotential is chosen so that the form of the Schwinger- 
Dyson equations will not depend on the choice of G (but the free energy will). 



^ The i in front of / gives the right sign for the kinetic term of /. 

^ This SU{N) becomes the color symmetry once we embed this theory into BFSS. 
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After integrating out the fermionic coordiantes, the action iq^ 

5 = / rfr Tr I \rr + \rJi + \rr 

' ' (2.4) 



where dot stands for Euchdean time derivative. The SUSY transformations are 

{Qa, 0] = -^^a, 
{Qa,i'p} = i^apf + i5ap(t), (2.5) 

[Qa,/] = -^ea/3V'/3- 

2.2 Schwinger-Dyson Equations 

The perturbative field theoretic approach to a (classically) massless theory in one dimension 
suffers from infrared divergences. In many examples, such infrared divergences are cured non- 
perturbatively. A framework that improves the ordinary perturbation theory and naturally 
resolves the IR divergences is the Schwinger-Dyson equations. While the usual formulation 
of Schwinger-Dyson equations are a set of recursive integral equations that express exact 
correlation functions or IPI vertices in terms of higher point vertices, such equations are 
often hard to solve due to the general momentum dependence in the exact vertices. In this 
paper, we will work with the equations that express the exact self-energies in terms of the 
integrals of exact propagators and tree-level vertices. We will make a truncation on the loop 
order of these equations, study the solutions to the truncated equations, and then discuss 
the validity of such truncations. 



Starting from the action (2.4), the Schwinger-Dyson equations at finite temperature are 
formulated as follows. We compactify Euclidean time on a circle of circumference /3, and 
expand the fields in their Kaluza-Klein modes along the Euclidean time circlq^ 






^° - y^ Z^ ^".-^ ' (2.6) 



rGZ+i 






2 

2'KinT I fi 



3 We let {T^} be a basis for SU{N) such that Tr iJ^T^) = 5^^, and write (j)" = (j)<'^T^. 
^ Unless otherwise noted, throughout this paper n, k, £ are integral and r, s are half-integral. 
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The action (2.4) written in terms of these Kaluza-Klein modes becomes 

—n^n ' cy / J n Ya,—rYa,r ' n / j J —nJ n 
r n 



^ = Tr^E'^'™^'- 



2^ /3 



+ 



IK 



2x/c^ 



n,k 



2Vc^ 



(2.7) 



Let us denote the exact propagators by 

\JnJm) = ^"" On,-m- 



(2i 



Note that {(f)f) is prohibited by the following Z2-symmetry (R-parity) of the action (2.4): 



—)■ 



(2.9) 



Working in 't Hooft units where the dimensionful 't Hooft coupling k'^N is set to 1, i.e. 
(n'^Ny^^jS — !■ /3, the "one-loop truncated" Schwinger-Dyson equations arq^ (see Figure [ij) 



1 

a: 



,27m ^2 2 Y^ 2 ^^ 



/3 



/3 



/3 



1 2nr 2 ^^ ^ 



(2.10) 



The terminology here requires some explanation. The equations (2.10) are the Schwinger- 



Dyson equations for the exact two-point functions, where the three- and higher-point func- 



tions are approximated by their bare values. We refer to (2.10) as the "one-loop truncated 



Schwinger-Dyson equations. The higher loop corrections will be discussed in Section 3.6 



This approach is somewhere in between a fully non-perturbative treatment and perturba- 
tion theory, but the results clearly go beyond perturbation theory since the nontrivial solu- 
tion for the self-energies cuts off the IR divergence that would invalidate the conventional 
perturbation theory. 

Our convention for the self-energies o"„, h^-, and ?7„, for the boson 0, fermion ■?/', and 

^ The sign for the fermion loop is compensated by the factor of {—if' from (-0^ ^V'fl s) = ^igrS°'^5ai35r.-s- 
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auxiliary field / respectively, is such that they are related to the exact propagators by 

1 



A, 



'" ( 2nn \ 2 



P 



+ <yn 



9: 



''- ^EL + h ' 
1 



(2.11) 



l + Vn' 



where a_„ = cr„, h-r = —K, and ?7_n = rjn- The Schwinger-Dyson equations (2.10), written 
in terms of the self-energies, are 

2 ^^ 1 2 ^^ 1 



CTr. 



(3 



E 



Vn=^Yl 



■ 2nk \2 
. /3 ) 



+ <^k 



[1 + r]n-k] 

1 



+ 7^E 



r \ o ~r lir 



27r(n— r) 



+ hn- 



(3 



[c-fy 



+ crk 



' 2n(n—k) \2 



+ Cr-n-k 



(2.12) 



(3 



E 



1 



[c-fy 



CTfe 



2n{r~k) 



+ hr 



2.3 J\f = 2 SUSY Ward Identities 

Although supersymmetry is broken at finite temperature, one expects supersymmetry to be 
"approximately" restored in the low temperature limit. This is subtle due to IR divergences, 
as the self-energies in the strict zero temperature limit are singular. Nonetheless, as shown 
later, the Ward identities we derive here assuming supersymmetry will hold approximately 
for the self-energies of the nonzero frequency modes of the fields in the low temperature 
limit. Such relations are useful in organizing the low temperature expansion of the solutions 



to the Schwinger-Dyson equations (2.12). 



Naively, at strictly zero temperature, where the momentum (in Euclidean time) is con- 



tinuous, we write the exact two-point functions (2.8) as 

{(l){p)(j){-p)) 



p^ + cr(p) ' 



(^a(p)^/3(-p)) 
(/(p)/(-p)) 



—i5. 



al3 



p + h{p) 
1 



(2.13) 



1 + ri{p) ■ 
If SUSY is not spontaneously broken, then 

= ({Q„,0(r)V^^(r')}) = -i{Mr)Mr'))+^5o.mrW)). 



(2.14) 
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Figure 1: The one-loop truncated Schwinger- Dyson equations with bare cubic couphng. 
The sohd, dashed, and double lines represent the boson propagator A„, the auxiliary field 
propagator e^, and the fermion propagator gr, respectively. 



where we used (0/) = 0. Similarly we have 

= ({Qa, /(r)^;3(r )}) = -ie^,{i^,{T)i^p{T')) + ze„^(/(r)/(r')). (2.15) 

In momentum space, (2.14) and ( |2.15[ ) read 

a{j)) = ph{p) = p^Tjij)). (2-16) 



At finite but low temperatures, if supersymmetry is only slightly broken, one may antic- 
ipate that there exists a continuous function s{p) such that 



(Tr, 



Vn 



,27m, 
,27rr 

■If' 



27m 



[I + 0(/3 
2Tcn 



+ 0{(5 



<-l^ 



<0 



(2.17) 



/3 



+ 0(/3 



<-l^ 



In particular, if we assume s(0) is finite, then one expects that ctq should be much smaller 
than (Jn^oi and r/o should be much larger than ?7„^o- This will indeed be the case. The 
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separation of zero modes and nonzero modes is the key to the scahng ansatz (3.1) that will 
allow us to solve the Schwinger- Dyson equations at low temperatures. 



3 The Low Temperature Limit 

3.1 A "Soft-Collinear" Approximation 

In the low temperature, or large /^ limit, we will demonstrate that the Schwinger-Dyson 



equations (2.12) admit solutions of the following leading order scaling in j3 (here the mode 



numbers n,r are assumed to be 0{1)) 

Note that the scaling of the zero modes ctq and r]Q are order (3~^^^ lower and higher, respec- 
tively, than their nonzero mode counterparts. This /3~^/^ turns out to be the appropriate 
expansion parameter for the Wess-Zumino quantum mechanics at low temperatures. 



One key property of the ansatz (3.1) is that at large (3 the boson zero mode propagator 



Ao ~ /3^/^ dominates the loop sums on the right-hand-side of the Schwinger-Dyson equations 



(2.12). Hence, to leading order in 1//3, the S-D equations reduce to simple algebraic expres- 
sions for the self-energies of the nonzero modes and for rjQ in terms of ctq- Once we have 
the leading order expressions, we can compute corrections in 1//3 order by order. Finally, 
substituting these expressions into the S-D equation for ctq, we find an equation involving 
ctq only. If a solution for o"o with the assumed scaling exists, then we know that the scaling 



ansatz (3.1 ) is indeed consistent. The consistency of the ansatz is highly nontrivial, however, 
and as we will see require delicate cancelations. By a slight abuse of terminology, we refer 
to our ansatz as a "soft coUinear" approximation, where we regard the nonzero modes as 
"hard" and zero modes as "soft" . 
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3.2 Nonzero Modes an, Tjni K and the Auxiliary Zero Mode tjq 

In the following we implicitly take n ^ 0. Let us reorganize the Schwinger-Dyson equations 



(2.12) for the nonzero modes an, Tjn, and hr and the auxiliary zero mode r/o in the form 



where 



_ 2 1 27ra 

fJao 1 + ?7„ (3 
_ 2 1 /3 

h -1 ^+c 



Vfe^O r 






Cr = -^ 7 ,^k9r-k, 

^ k^O 



According to the scaling ansatz (3.1) 



' ' r^'\ ^A.~/3-^ 



(3.2) 



(3.3) 



(3.4) 



/3ao l + Tjn (3 

H /27rn\9 I ~ /3 ) 7, Bn-^0^'^ 

Since the terms on the RHS involving A^, Bn, Cr, and D are sub leading in 1//3, the leading 



order Schwinger-Dyson equations (3.2) are simple algebraic equations as advertised earlier. 
Notice that because we expanded in 1//3 while taking n,r ^^ C'(l), the loop sums where n 
and r range over — 00 to cxd become divergent and will require regularization. 
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The fermion self-energy hr and the product ant]n can be directly read off 



sign(r)^^ + 0(/3-'/'), (3.5) 

anVn = ^ + Oir'^'), (3.6) 

PCTO 

where sign(r) conies from the odd parity of hr. To find an and ?7„ individually, let 

2?™ 



27m 

where 



"^nt-L ~r Xfi) , 



(3.7) 



s„^sign(n)^/— + 0(/3-2/5). (3.^ 



If the approximate SUSY Ward identities (2.16) are to hold in the low temperature limit, 
then we expect x„ to scale with a negative power in /3. In Appendix [D] we solve for Xn 
directly from the Schwinger-Dyson equations, and find 

Xn = i^)^w^, + oir'^') - /3-^/^ (3.9) 

z Z7i\n\ 



in accordance with (2.16). 



3.3 Boson Zero Mode ctq 

Let us examine the Schwinger-Dyson equation for the boson zero mode ctq 

If one substitutes the leading order expressions for the other self-energies into the above 
equation, then each term becomes order (3~'^^^, and one finds that the two sides exactly 
cancel. It turns out that the next order /3~^ terms also cancel, and therefore the solution to 
(Jo is determined at order f3~^'^^^. This requires expanding s„ and h^ to order (3~^ and tjq to 
order (3^^^ (each two orders of (3^^^^ down from the leading piece). 

For the auxiliary zero mode rjo, the relevant correction is given straightforwardly by D 



in (3.2), and we find 



*-^4E45^il(ri^ + °(r-)-^^^ + °(A ,3.u; 
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For Sn and hr, let us write 



hr = sign(r)y^ + /.W + /i(2) + 0(/3-8/5), 
where Sn , hr are of order /3~^/^ and Sn , hr are of order (3~^. Next, define 



(3.12) 



9 



l^^o ^ _ v^ sign (£) sign (A; 



(3.13) 



e^o,k 



r ) 



where k can take both integral and half-integral values. The computation for Sn , Sn , h 
and hr can be found in Appendix [EJ The results are 



/,(!) = ^Cr + 0{P~'] 



Air 



^^ = -^ + sign(n)/^C^-f^ 






+ 



" E (-)'-g"(« - ^)c'^/- i E 



47r2n 



fe^0,2n 



E 

2 



sign(n — k)Ck + sign(A;)C, 



n—k 



k 



^ — sign(A;)sign(n — k)xk 

2tt 



k=/=0,n 



k 



+ 



h^ 



Tvr 
'J 



+ sign(r)#-a2 



0(/3-8/5), 
^ (7^ Y^ sign(r — k)Ck + sign(A;)Cr 



in 



2 r 



Att- 



E 

k^O 



k 



g^ \--^ sign(fc)sign(r — k)xk 



2n 



E 

k^O 



k 



0{(3 



-8/5N 



The relevant contribution to the Schwinger- Dyson equation for o"o is 



n^O r 



-'- -^M 1 + ^ + ^ ) + o{P'''"^ 



Using (3.11) and (3.15), the equation for ctq becomes 



(3.14) 



(3.15) 



^0 = ao + :^/o"p-"' - YqP'4 + 0{p-^''% 



from which we obtain 



cro = 2(0'/3-i + 0(/3-^). 



(3.16) 
(3.17) 



To summarize, we solved the Schwinger-Dyson equations and obtained the self-energies 



at low temperatures. Consistency of the ansatz (3.1) is thus verified. Note that although we 



solved Sn to second order in (3.14), x„ is only computed to zeroth order in (3.9), so according 



to (3.7) we only have first order expressions for an and ?7„. Fortunately, for the purpose of 



solving (To and later on computing the free energy, we only need Sn to second order but not 
an and rjn separately. 

3.4 Continuum Limit 

Before computing the free energy, let us note the following property of our solution of self- 
energies. In terms of the low temperature expansion parameter g, the boson self-energies 
are 



o-Q 



an 



27r|r7,| 



1 



27i\n\ 



2n 



(3.18) 



0(/3 



-2\ 



where the first equation is the definition of g. If we write p = 27rn//3, then the self-enengy 
for the nonzero modes can be written as 



a{p) 



M 

9 

M 

9 



1 + 



3^ 



1 + 



/3 



\n\ 



2E 



fc=l 



2nk 



+ 0{/3- 



(3.19) 



1 + 



2^ 



9' 



P\p\- 



/3\p\ TT ' 27r 



+ 0{g' 



where HN{x) is the analytic continuation of the harmonic function. Since HN{x) is regular 
at a; = 0, taking p — )■ gives 

29' 



a{0) 



a 



+ 0{g'). 



(3.20) 



Curiously, (t(0) coincides with ctq. In other words, the zero mode self-energy a^ could be 
regarded as the p — )■ limit of the self-energy function a{p) of a continuous momentum p. 
This fact may be useful in applying the Schwinger-Dyson equation to Wess-Zumino quantum 
mechanics coupled to gauge fields, where a continuous distribution of the eigenvalues of the 
holonomy matrix along the thermal circle is involved. 

3.5 Free Energy 

Having obtained the propagators from the one-loop truncated Schwinger-Dyson equation, 
we can now compute the free energy in the mean-field approximation [23]. Let us define a 
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quadratic Fourier space "trial action" 5*0, whose propagators are solutions to the one- loop 
truncated Schwinger-Dyson equations: 






a 
nJ —n 



2gr 



2e. 



(3.2i: 



Next we write the true action as S* = 5*2 + S'3, where S'2 and ^3 consist of the quadratic and 
cubic terms, respectively. Then an approximate free energy is given by 



/3F = /3Fo + (^2-^o)o-2(^3')o. 



(3.22) 



The subscript here indicates that the expectation value is taken with respect the trial 
action, i.e., 

{■■■)o = J D<j)D^Dfe-^'^--- . (3.23) 

In the planar limit, and dropping an overall factor of N'^Nf, we have 
/3^o = -J$^lnA„-^5^1ne„ + ^J]ln^,^ 



(S'2 - S'o)o - 7,/ , 



2 

] 
2/3 



/27m\ , 



+ ^E(^"-1)-E(^9--i) 



(3.24) 



1 1 1 

9('S'3)o = W^ 2^ AmAfcem+fc + -^ / ^grQs^r 



m,k r,s 

In Appendix [F], we compute these contributions to the free energy using solutions computed 



in Section 3.3 but without plugging in the explicit value of ctq in (3.17). The results are 



Vn^O r / 



' \n>0 r>0 



192 



(/3ao)3 + 0(/3-9/^) 



(52-So)o = l-^(^)^/^ + 0(r^/^), 

\n^0 r / 



(3.25) 



/3 



E^-Em + t^(/5^o)' + ^(/5"'^')- 



,n>0 r>0 
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We see that the potentially divergent sums cancel between jSFq and -~i^{S^)o, and thus 

PF = const + l(/3ao)^ - ^(^)'/' + 0{r'/% (3.26) 

In |23], the one-loop truncated Schwinger-Dyson equations are obtained by extremizing 



the mean-field free energy (3.22) with respect to the propagators in the trial action (3.21) 



(which is why it was called the trial action in the first place). As a consistency check, we 



should find that our solution for o"o (3.17) minimizes the free energy (3.26). This is indeed 



the case as one can easily verify. Replacing ctq by its explicit solution (3.17), we get 



PF = const - ^(|)^/^Ar2Ar^/3-6/5 + 0{(3-''^'), (3.27) 

where we have restored the overall N'^Nf factor. 

To conclude, we solved the one-loop truncated Schwinger-Dyson equations in the large /3 
limit, and found that the non-constant part of /3F has (3~^^^ scaling. We have also confirmed 
it numerically, including the coefficient in front of /3~^/^. For comparison, the numerical 
analysis of [23] by fitting up to /3 ~ 28 gave the scaling ~ /3~^-°^. In fact, a careful inspection 
of various terms in the ctq S-D equation shows that the correct /3~^'^ low temperature scaling 
is visible only when /3 is greater than ~ 100 (in the normalization convention of [2^). 



3.6 Higher Loop Corrections 

A priori, there seems to be no reason to expect the one-loop truncation of the Schwinger- 
Dyson equations to capture the correct low temperature physics, since the dimensionless 't 
Hooft coupling is strong in the low temperature limit. Naively, one may expect two and 
higher loop contributions to the S-D equations to overwhelm the one-loop contribution and 
completely destroy the scaling behavior we found. However, due to the peculiar scaling 
behavior of the propagators, the effective expansion parameter in the low temperature limit 
is not the 't Hooft coupling. If we proceed with our ansatz for the one-loop self-energies, the 
question is whether the correction to the S-D equation for o"o respects the scaling ansatz. 

Let us look at the two-loop truncated Schwinger-Dyson equations, which are shown 
diagrammatically in Figure [2] As it turns out, the two-loop diagrams contribute an additional 
terno 

^ (^ - ^.\ f^l ''' In A ^ o{p-- In /3) (3.28) 



The log comes from cutting off the harmonic sum at ^/Jija^^. 
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Figure 2: The two-loop truncated Schwinger-Dyson equations with bare cubic couphng. 
The sohd, dashed, and double lines represent the boson propagator A„, the auxiliary field 
propagator e^, and the fermion propagator g^, respectively. 



to the RHS of the Schwinger-Dyson equation for ctq (3.16 ). This is, up to a logarithmic factor 



the same as the next-to-leading order term in the one-loop contribution. However, as seen 
earlier, the scaling solution to the one-loop S-D equations requires complete cancelation at 
next-to-leading order, and ctq ~ (3~'^^^ would be determined only at next-to-next-to-leading 
order. This is now spoiled by the two-loop contribution. This means that either the ansatz 



(3.1) is incorrect, or that x„ is at least of order 0{l) which invalidates the approximate 



SUSY Ward identities (2.17). So it appears that the two-loop truncated S-D equation of the 
M = 2 Wess-Zumino quantum mechanics does not give rise to a simple scaling solution in 
the low temperature limit. 

Inclusion of three- and higher-loop diagrams would give corrections that even overwhelm 
the one- loop leading order terms in the low temperature limit, and hence completely invali- 



dates the ansatz (3.1). Figure k5^ shows two diagrams appearing in the pertubative expansion 



of the exact three-point vertex that are larger than their tree-level counterparts. Diagram (a) 
does not exist at tree-level, and Diagram (b) is of order 



P'^'I^A^ftn ~ r^'\P'/''fP'^"' ~ P-^'^P"\ 



(3.29) 
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which is (3'^/^ larger than tree-level (/3~^/^). It is possible that a scaling behavior of the type 
we found still holds in the full theory at large A^ jj but if so, an alternative truncation of the 
Schwinger-Dyson equation would be needed to reveal the correct low temperature scaling 
beyond the one-loop result. 





(a) (b) 



Figure 3: Diagrams that appear in the pertubative expansion of the exact three-point 
vertex. 



4 A/^ = 4 Wess-Zumino Matrix Quantum Mechanics 

The A/" = 4 Wess-Zumino matrix quantum mechanics, whose one-loop truncated Schwinger- 
Dyson equations take an identical form as those of the M = 2 theory, on the other hand 
does not receive two-loop contributions, nor planar three-loop contributions. Moreover, the 
perturbative expansion of the exact three-point vertex only contains diagrams that are at 
most of the same order as the leading order at one- loop. It is therefore conceivable that 
one-loop truncated S-D equations serve as a better approximation of the Af = 4 model. An 
A/" = 4 Wess-Zumino matrix quantum mechanics can also be obtained as a truncation (or 
deformation) of the BFSS matrix model. In the N' = 4 language, the BFSS matrix model 
consists of one vector multiplet and three chiral multipletsjj with an SU{3) invariant cubic 
potential. 

^We will comment on the relevance of the large N, i.e. planar, limit in the next section. 
* Since ID A/" = 4 is the dimensional reduction of 4D Af = 1, wc name the supermultiplets using the 4D 
language. 
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4.1 Action, Schwinger-Dyson Equations, and Free Energy 

Consider an A/" = 4 Wess-Zumino matrix quantum mechanics with a flavor symmetry group 
G. Suppose the chiral multiplets are in an A^/-dimensional representation of G, governed 



by a cubic superpotential of the same form as in explained at the beginning of Section 2.1 
(except that the real M = 2 matter multiplets are now replaced by the holomorphic A/" = 4 
chiral multiplets). For the A^ = 4 WZ model coming from the matter sector of BFSS, 
G = SU{3) and Nf = 3, c = 2. Our conventions for the ID A/" = 4 superspace are defined 
in Appendix O We have Nf chiral superfields 



<!>- = r{y) + V2ra{y)9'^ + zny)9' 



(4.1) 



in the adjoint representation of SU{N) and in the fundamental representation of the flavor 
symmetry G. The superpotential is given by 



W 



6v^ 



fte'^^'^Tr ($'^[$^$"]) 



and the action is 



S = / dr Tr <^ 0" + ^""C + FT 



+ I ^fn',<P'] + ^e"^e"V[V^',^^] +h.c. 



(4.2) 



(4.3) 



Let us write the finite temperature propagators in momentum (frequency) space as 



'^:^l 



) = A„5" 6n. 



\JnJml — *="" '-'n-m 



(4.4) 



The self-energies are defined as in (2.11). The A/" = 4 SUSY Ward identities, given in 



(C.24), take an identical form as the Af = 2 identities (2.16). We have normalized the 



coupling constant n in (4.3) such that the one-loop truncated Schwinger-Dyson equations 



are identical to the Af = 2 ones (2.12), whose solution is given in Section 3.3 



The free energy can also be computed to be 

5 ,7C , 



(3F = const - -{-f/^N^NfP-^/^ + 0(/3-^/^). 



(4.5) 



which is twice that of the M = 2 free energy (3.27) because the fields are now complex. To 
restore the dependence on k^, one replaces (3 by the dimensionless combination (3{k'^NY/^. 
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4.2 Higher Loops 

A key difference between the A/" = 2 and A/" = 4 Wess-Zumino quantum mechanics is the 
holomorphicity of the superpotential. In the Feynman diagrams of the A/" = 4 model, propa- 
gators carry arrows since the fields are complex, and vertices have either all three propagators 
pointing inwards or all pointing outwards; in other words, the Feynman diagrams are bipar- 
tite graphs. It is not hard to see that no two-loop or planar three-loop diagrams can appear 
on the RHS of the Schwinger-Dyson equations for the self-energies. Thus without taking the 
planar limit, the first nontrivial corrections come from non-planar three-loop diagrams, and 
in the planar limit, the first nontrivial corrections come from four-loop diagrams. Examples 
of these non-planar three-loop and planar four-loop diagrams are shown in Figure |4} 



.-r 



>- 



(b) 



Figure 4: Corrections to the auxiliary self-energy rjn from (a) a non-planar three-loop diagram 
and (b) a planar four-loop diagram. 



Although we have not explicitly evaluated the four-loop contributions and its effect on 
the scaling solution in the low temperature limit, a naive counting of powers of propagators 
of zero modes versus nonzero modes based on the one-loop scaling ansatz suggests that the 
higher loop contributions are of the same order as the one-loop contribution. While such 
contributions are expected to correct the coefficient of the self-energies and the free energy, 
it is conceivable that the scaling behavior found in the solution to the one-loop truncated 
equations continue to hold when all-loop contributions are included, at the planar level. The 
large N limit is potentially important here, in order for the loop expansion in the Schwinger- 
Dyson equation to converge. At L-loop order, there are roughly L\ 0{C^) diagrams that 
contribute to the Schwinger-Dyson equations at leading order in 1//3, but only 0{C'^) planar 
diagrams that do so, where C and C are some constants [29] . For example. Figure IS] shows 



an L-loop correction to the fermion self-energy hr. Assuming the scalings in (3.1 ), the scaling 
of this diagram is 



r'^g^-'A^g^ ~ /3-^/3 



Ln(-L+l)/5 



/3^^/5/3- 



L/5 



/3 



1/5 



(4.6) 
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which indeed is the leading order in the Schwinger-Dyson equation for hr. There are two such 
planar diagrams, coming from the two (alternating) assignments of the R-symmetry index 
q; = 1, 2 at the vertices in the fermion loop. On the other hand, there are 2(L/2)!(L/2 — 1)! 
non-planar diagrams, where the extra (L/2)!(L/2 — 1)! comes from inequivalent permuta- 
tions among the boson propagators. While summing over all non-planar diagrams gives an 
asymptotic series, restricting to planar diagrams may give rise to a convergent sum. 




Figure 5: A leading order correction to the fermion self-energy hr from a planar L-loop 
diagram {L has to be even). There are L boson zero mode propagators (solid line). 

We conjecture that the /3~®/^-scaling of the free energy is preserved to all loop orders for 
the A/" = 4 Wess-Zumino quantum mechanics in the planar limit, with only the coefficient 
to be corrected. Further, it is possible that such scaling behavior holds when perturbative 
1/N corrections are included as well, and is only violated by non-perturbative effects in 
1/A^. One reason to anticipate this is the existence of exponentially long lived metastable 
states in BFSS matrix quantum mechanics, which dominate the thermodyamics as predicted 
by the gravity dual. As explained in the introduction section, the scaling behavior cannot 
possibly hold for the exact free energy in the general Wess-Zumino quantum mechanics at 
finite N: depending on the precise matter content and the form of the superpotential, the 
model could be either gapped, in which case the free energy is exponentially suppressed 
in the very low temperature limit (for fixed N), or infinite, when there is a moduli space 
(fiat directions of the superpotential). Nonetheless, in the infinite A^ limit, we expect that 
there is at least a phase of the theory in which the free energy exhibits the nontrivial power 
scaling in the temperature, as we found from the one-loop Schwinger-Dyson equation. In 
A/" = 4 theories where the superpotential is such that the spectrum is gapped at finite A^, we 
expect a near-continuum is developed in the large A^ limit, and if we first take A^ — > oo, and 
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then take the temperature to be small in 't Hooft units, the free energy should exhibit the 
scaling behavior. In A/" = 4 theories where the superpotential has flat directions, there are 
alternative, and singular, solutions to the Schwinger-Dyson equation, where the self-energies 
of the bosonic fields corresponding to the flat directions are set to zero, while the self-energies 
of the fields along the non-fiat directions diverge. Such singular solutions may be regularized 
by introducing an IR regulator that cuts off the volume divergence in the target space. 
In such theories, at infinite A^, we expect our scaling solution to describe an "unbroken" 
phase (the flavor group is unbroken), as opposed to the "broken" phase describing the flat 
directions of the potential. The different phases correspond to different solutions to the 
Schwinger-Dyson equations. 



5 Towards the Low Temperature Expansion of BFSS 

In the previous sections we have studied the low temperature expansion of the M = 2 
and A/" = 4 Wess-Zumino quantum mechanics. To extend our results for the Wess-Zumino 
quantum mechanics to BFSS matrix theory, we need to couple the vector multiplet to it 
and perform the appropriate supersymmetric gauge fixing, which also introduces the ghost 
multiplet. In this section we discuss some preliminary results toward understanding the 
planar free energy of BFSS. 

5.1 A/" = 2 Gauge-Fixed BFSS 

The BFSS action can be decomposed into one M = 2 vector multiplet F^ plus seven matter 
multiplets $", a = 1, ■ ■ ■ , 7 in the fundamental representation of the G^ flavor symmetry. 



Svec + Smaller =^ / dTcfOTl | " T ( ^(^a/S j VaJ^'V/sJ^'j 



(5.1) 



where J-"* is the field-strength superfield for the real connection superfield Fq,. See Appendix 
[B]for our notations. 

The Af = 2 gauge-fixing condition is 

e^pD^r^ = 0, (5.2) 

If we expand the real connection superfield Fq, as 

Ta = Xa + iAoOa + X'^^^Ofs + dea^Op + 2eaisXisO^. (5.3) 
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(5.2) becomes 



Ao = 0, d = 0, Xa = -^Xa. (5.4) 
The Af = 2 gauge-fixing introduces a ghost multiplet 

C = a + PJc + ilO'^ (5.5) 

with a and 7 complex fermonic fields and /3q a complex bosonic field. The ghost action is 

Sghost = fdrdO' Ti Ue^pD^CVf}C\. (5.6) 



(5.1) together with (5.6) gives the complete BFSS action in the M = 2 gauge (5.2). 



5.2 The M =2 Ghost Determinant 

In the usual Wess-Zumino gauge-fixing, which breaks all supersymmetries, the eigenvalues of 
the zero mode of Aq have the 2tt periodicity coming from the large gauge transformation. Or 
equivalently, the eigenvalues of the holonomy matrix U = exp (^ drAoj are gauge invariant. 
However, there is no such periodicity for the Af = 2 gauge-fixing. To see this, let us note 
the coupling between X\ f3a, and the gauge field Aq in the BFSS action. 



il\2 



S,ec3{drX' + t[Ao,X']) 
Sghost 3 Pa (Pa + ^[^o. Pa]] + \liMX\ Ppi 

where 7^ = cri,7^ = (T3. Note the curious 1/2 factor in the coupling between the /3-ghost 
and the gauge field. Suppose we shift one of the eigenvalues of Aq by 47r, then this can be 
undone by shifting the momentum modes of X* and Pa by 2 and 1, respectively. However, 
the coupling between X* and Pa is not invariant under this shift, so the 1/2 factor in the 
gauge coupling of Pa completely destroys the periodicity of the eigenvalues of Aq. 

In fact, the J\f = 2 gauge condition (5.2), ( |5.4 ) leaves no residual large gauge transfor- 



mations (except for the constant gauge rotation). First consider the zero temperature case. 
Recall the infinitesimal M = 2 gauge transformations form Appendix A. 2.1, generated by 



A = f2 + icUaOa + W9 . The ones that preserve the gauge condition (5.2) obey 



d f ■ i 

0= — (n + t[Ao,n] + -{xa,(^a} 

= -W-^eap{Xa,oJp}, (5.8) 

i 1 i d 1 

= UJa + ^[Ao,UJa] + -eap[Xp,W] - [A„, fi] + 2d^[Xa,^ " "7^/3 [^\ ^/s] • 
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The last equation can also be written as 



i 1 i ■ 1 ■ ■ 



(5.9) 



For generic Xa, these conditions are satisfied only if Ua = 0, W = 0, and ^2 = 0. Therefore, 
the residual unfixed gauge transformations are constant bosonic gauge rotations. This is 
in contrast with the non-supersymmetric gauge fixing (Wess-Zumino gauge together with 
Aq = 0) where Q can be linear in time r, provided that it also commutes with Aq. In 
the finite temperature case, the finite gauge transformations must be such that the bosonic 
components of e^ is periodic on the Euclidean time circle while the fermionic components 
are anti-periodic. While in the case of non-supersymmetric gauge fixing, there are residual 
large gauge transformations of the form e*'''^ with [Ao,C] = and e*^^ = 1, in the case of 
M = 2 gauge fixing the only residual gauge transformations are e^ with constant bosonic Vt. 
Therefore, in the latter case, the eigenvalues of Aq are not periodically identified, and in the 
path integral we must integrate over all range of Aq. 

If we set Xa to zero, the M = 2 ghost determinant as a functional of (constant) Aq and 

X^ is 

det\dr + i[AQ,-]) 

det {{d, + |[Ao, ■])5^p + \l'^^AX\ ■]) ■ ^ ^ 

As already mentioned, if we further restrict to X* = 0, this ghost determinant is invariant 
under the shift of eigenvalues of ^Aq by 47r, but not by 27r. The path integral measure in 
Aq can be written in terms of the integration over its eigenvalues ctj (j = 1, 2, ■ ■ ■ , A^) as 



N 



JJrffti JJ 



1=1 



i<3 



sin^(Q;ij/2) 
cos4(aij/4)' 



(5.11) 



In the large A^ limit, in terms of the eigenvalue distribution function p{a) (whose integral is 
normalized to 1), the effective potential for a due to the ghost determinant is (when Xa and 
X* are restricted to zero) N"^ times 



1 



dada p{a)p{a) 



— In sin 
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-i)V 



n/2 



n=l 



n=l 



n 



n=0 ^ 



(5.12) 



where 



Pn 



dap{a) cos{na). 



(5.13) 



Note that the coupling of matter multiplets to the vector multiplet still respects the shift of 
f3AQ by 2tt, and consequently, the effective action due to integrating out matter multiplets 
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should depend only on p„ for integer n and not p„+i. To minimize the effective action, one 
would set p„_|_i to zero for all integer n, and effectively the eigenvalues no longer repel one 
another. This is rather surprising. Let us note cautiously, however, that we have essen- 
tially ignored the scalar X* and fermion Xa in the vector multiplet, in the above discussion. 
Integrating them out could have important consequences on the effective action of Aq. 

5.3 The Free Energy 

To obtain the free energy for BFSS, one needs to incorporate the vector sector in a con- 
sistent manner. For example, one could try to derive the effective action for the vector 
multiplet by integrating out the matter multiplets, using the solutions to the one-loop trun- 
cated Schwinger-Dyson equations of the Wess-Zumino quantum mechanics, around a certain 
background of Aq. It turns out that the Schwinger-Dyson equations are not particularly use- 
ful in solving the vector sector with the effective potential l^(Ao) generated by the matter 
sector. The difficulty lies in that, the soft collinear approximation doesn't hold for the vector 
and ghost multiplet. Nonetheless, we can infer some properties of V{Aq) from our results. 

Assuming that the matter sector has the same scaling behavior as in the M = 2 Wess- 
Zumino matrix quantum mechanics, integrating out the matter sector to leading order in 
P~^ generates an effective mass term for the vector multiplet 



my 



j dr Tr {Al + {X^f + x.Xa) , < ^ ^ ~ ^'^'- (5-14) 



In the more familiar case of non-supersymmetric gauge fixing, while typically the effective 
potential for Aq drives the eigenvalues of Aq toward the origin, there is a competing repelling 
effect between the eigenvalues due to the ghost determinant (or the measure in the integration 
over the eigenvalues). Naively, the effective mass my would be too small to overcome the 
repelling of the eigenvalues. However, as seen in the previous subsection, in the J\f = 2 
gauge fixed path integral, the eigenvalues of Aq in fact do not repel, as far as the coupling 
to matter multiplets is concerned, as the latter only knows about the eigenvalues of Aq 
up to shifts by 27r//3. It then seems plausible that the free energy of the vector multiplet 
coupled to matter multiplets is dominated by the contribution from the matter multiplets 
alone, with the vector multiplet turned off. If this is the case, the free energy of BFSS 
matrix quantum mechanics would have the same low temperature scaling as the truncated 
Wess-Zumino quantum mechanics, and we would have (3F ~ —CT^^^. We do not know the 
source of the discrepancy between this 6/5 scaling exponent and the 9/5 as predicted from 
the gravity dual. Clearly, a better understanding of the dynamics of the gauge multiplet, 
through perhaps a better set of truncated Schwinger-Dyson equations that do capture the 
correct low temperature behavior, is needed. 
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Let us also note that, unlike the typical large A^ gauged matrix models with adjoint 
matter, where one expects a Gross- Witten-Wadia phase transition [301 W^ going from high 
temperatures to low temperatures, such a transition should be absent in the planar limit of 
BFSS matrix quantum mechanics, because there is no Hawking-Page transition of the DO 
black hole in the gravity dual, for temperatures of order 1 in units of the 't Hooft coupling. 

Previously, the authors of [221 [211 125] studied the Schwinger-Dyson equations of this 
Af = 2 gauge- fixed BFSS model (presented in Appendix B.4), and produced numerical 
results for the free energy. A /3F ~ /J"^'"^ scaling was found there, which appeared to be 
close to the /3~^-^ predicted by the dual black hole. However, this scaling was obtained by 
fitting the numerical results within the inverse temperature range 1 < /3 < 4. If one modifies 
the fitting range to, say, 0.25 < /3 < 4, the result already changes drastically. In addition, 
the numerics break down at /3 ~ 6, essentially due to the ghost multiplet. In fact, the 



Schwinger-Dyson equations for just the pure gauge system (B.57), (B.56) do not admit a 



real solution due to the negative signs in the equations for 11^ and S^. When /3 > 6, the 
terms from coupling to matter become smaller than the pure gauge terms, and it appears 
that no real solution exists in this regime. Further, it is unclear to us why [211 EH] treated 
the gauge field as a holonomy matrix model after gauge-fixing and including the ghosts in 
the S-D equations: such a treatment would seem to double-count the ghost determinant. 
Though, this double-counting turns out to have little effect on the numerical result for the 
free energy computed from the S-D equations in the temperature range considered in 



5.4 J\f = 4 Gauge-Fixing 

As mentioned before, the Af = 4 Wess-Zumino matrix quantum mechanics with SU{3) fiavor 
can be regarded as a deformation of BFSS. In A/" = 4 language, the BFSS action can be 
written as 



£ 



BFSS 



Tr 



1 



4^ 



2 
YM 



z'^^WJV. 



+ h.c. 



+ 



<J,ag-2y^a 



D 



6^2 



^yAfe"'''^<l>"[<l>^<l>" 



h.c. 



(5.15) 



where \^ is a real superfield representing the vector multiplet and W is the field strength. 
The advantage of writing the BFSS action in A/" = 4 language is that we can impose a 
manifestly A/" = 4 supersymmetric gauge fixing condition 



e^^V^V^V = 0. 



(5.16) 
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where V is the real superfield representing the vector muhiplet. Pq, is defined as 



'^a = -l^,-d^S,^dr (5.17) 



Expanding the real superfield V{t, 9, 0) as 

V(r, e, 6) = R{t) + Oxir) + ^x(r) + OOMir) + MM^t) 



(5.1^ 



then (5.16) implies 

M = M = 0, 



K = -^^apX^, K = T^SpaX^ (5.19) 



D = Ao- -R. 

" 2 



Note that since the lowest component field R{t) is not zero in the A/" = 4 supersymmetric 



gauge- fixing condition (5.16), there are infinitely many terms involving R{t) in the action. 

For now, we wish to decouple the vector multiplet and only study the matter (chiral) 
multiplets described by the Af = 4 Wess-Zumino quantum mechanics. This can be done by 
deforming the BFSS action with a superpotential SW = -=^Ke"*'^Tr ($"[$'', •l)'^]) and taking 
i^/qym -^ oo. Note that this deformation is in the same supermultiplet as Tr ($"$"), which 
is presumably dual to a stringy mode in the bulk. Hence the gravity dual of the deformed 
theory has a very stringy description. 

5.5 Phases of BFSS 

There can be different solutions to the one-loop truncated Schwinger-Dyson equations. In 
Section |4| we have only considered the solution to the A/" = 4 matter sector WZ model that 
respects the SU{3) fiavor symmetry. To obtain a solution where the SU{3) flavor symmetry 
is spontaneously broken, we can take the zero mode self-energy ctq of, say, $^ to be very 
small, which corresponds to $^ being massless and can acquire a large vacuum expectation 
value. This in turn generates large masses for $^ and $^. We can then integrate out $^ 
and $^ and obtain some effective action for $^. In the effective action there are scattering 
states described by wave packets at large $^, that becomes just that of the free theory of one 
chiral superfleld $^. The free energy is then given by A^^ times the volume divergence from 
the scattering states, which is certainly different from the /3-scaling we found in Section |4j 
As already discussed, we expect these different solutions to the planar Schwinger-Dyson 
equation to describe different phases of the theory at inflnite N; the tunneling from one 
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phase to another, which resembles the Hawking radiation of DO black holes in the gravity 
dual of BFSS matrix theory, is expected to be exponentially suppressed in the large N limit. 

When the matter WZ model is coupled to the gauge multiplets, the ~ A^^ flat directions 
are reduced to ~ A^ fiat directions, due to gauge symmetry. From the perspective of the 
planar Schwinger-Dyson equation, there is no longer the singular solution that breaks the 
SU{3) flavor symmetry, as the loops of gauge multiplets contribute to the self-energy of all 
matter chiral multiplets. Given that the planar free energy is supposed to capture metastable 
microstates of the DO black hole, an intriguing question is whether one can study quantita- 
tively the tunneling amplitudes which are non-perturbative in 1/A^, in the framework of the 
Schwinger-Dyson equations. This is left for the future. 

6 The High Temperature Limit 

Compared to the low temperature analysis of BFSS in the previous sections, the Schwinger- 
Dyson equations in the high temperature limit are much simpler to deal with, because the 
nonzero modes are weakly coupled (the zero modes remain strongly coupled). Prior work in 
this aspect includes [32], where the high temperature behavior of BFSS was studied by first 
dimensionally reducing to IKKT matrix model, and then utilizing Monte Carlo method to 
extract the two point functions. Here, we shall employ the one-loop truncated Schwinger- 
Dyson equations to perform the integrating out procedure. 

While the low temperature limit of BFSS matrix model describes a black hole in the bulk, 
its high temperature limit corresponds to the stringy regime in the dual theory. Therefore 
the high temperature analysis should tells us about the stringy black hole in the bulk. Note 
that while the horizon size of such stringy black holes are small in units of the string length, 
they are large with respect to the Planck length determined by the local string coupling 
(which goes to zero at large radii). 

For the following analysis, we will work with the M = 2 gauge fixing, solve the one- 
loop truncated Schwinger-Dyson equations in the /3 ^ 1 limit pertubatively, while keeping 
the next to leading corrections to the self-energies of both zero and nonzero modes (the 
appropriate expansion parameter in the high temperature turns out to be /3^/^), and obtain 
the free energy up to 0{(5^). Note that to obtain the next to leading correction to the free 
energy, we only need to know the /3^/^ corrections for the zero mode self-energies. We shall 
also calculate the "size" of the bound state (-R^)^, and compare with the Monte Carlo result 
from 
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6.1 Schwinger-Dyson Equations 

At leading order, the Schwinger-Dyson equations reduce to those involving only the zero 
modes, and those which determine the nonzero mode self-energies from the zero mode self- 
energies. The equations for zero modes at leading order are 

nM = ^ + A + 12 



^y _ 2 14 2 



u 2 7 

+ 



The detailed calculation is presented in Appendix G 



(6.1) 



2 m^ /3n^' 
whose solution is given by 

with 

b = c = 4.78511, a = ^^—77 = 3.54504, d = -^ = 0.477429. (6.3) 

b — 4/0 a'' 

Then one can readily obtain the leading self-energies for the nonzero modes 



Boson : H^ ^ (f - %)P''/\ K ^ (b - l)P~"\ Hf ^ ^J'l\ 

Fermion: E^ ^ ^^(6 - i)r'/2, Sf ^ ^f (6 - |;)r^/2, (6.4) 

Ghost: n„«=.-ir^^ i:^^i,ip--l\ 2„^-Te^/3^. 

Proceeding to next order in /3^/^, we find for zero modes, 

nj,'^"' = -2.52/3, ni')'' = -0.32/3, 'e!^^'' = Qmi3^/\ //« = 5.53/3, (6.5) 

while for nonzero modes, we have 

r 'n("" = 5.24183 + '^^, r 'n™-- = 4.86502 + ^H?1H. (6.6) 

^Z ^Z 
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6.2 The Free Energy 

Now that we have the self-energies, we can immediately obtain the free energy in the high 
temperature limit. To leading order, we have 

/3F ~ const + 6 log/3. (6.7) 

This agrees with the numerical result in [25] . 

Taking into account the nonzero modes, and also the subleading correction to the zero 
mode self-energies, we compute the free energy to the subheading order 0(/3^/^) in Ap- 



pendix G 



/3F = const + 6 log/3 - 3.89872/3^/2 ^ ^^^^3^ ^g_g^ 

6.3 Size of the Wave Function 

From the knowledge about the self-energies of the scalar fields, we may calculate the observ- 
able that measures the size of the wave functionP 



{R')p ^ ^ (Tr X^X^)^ + -^ (Tr r^)^, (6-9) 



where 

^(TrX^X% = 2r'J]a^, 

n 

n 

which to the leading order in high temperature expansion become 



(6.10) 



^(TrA-A-),.^ = ! 



1 „ .. ... 7 7._,„ ,6 2, 



^.(T.rr).-^ = -r-/^ = (^-^)r-/^ 



(6.ii; 



Hence we have 



{R'), ^ l(3-'%'yMN)'^ = 2.392gyMN'^'T'/', (6.12) 

where we have restored the 't Hooft coupling gy^^^ which was set to 1 previously, by 
dimensional analysiq^ In the dual supergravity geometry, the IIA DO black hole has radius 



^ Note that our normalization for the fields differs from 32 by a factor of y/N. 
^'^Recall that X^ has mass dimension 1 and gyM ^^^ mass dimension 3. 
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r ~ iV^/^T^/^ in the gvM = 1 unit, which agrees with the spread of the wave function ^J~{B?)p 
in the 't Hooft region T ~ N^l'^. 

In Appendix G we continue to compute {R'^)i3 to the subleading order 0{P), 



(6.13) 
where 



{R')^ = u,P^ 


-1/2 


+ U2l3 + OiP'/^) 


b 

- = 2.392, 


UJ2 


2/3 12 



coi = - = 2.392, u;2 = j^-- = 2.180. (6.14) 

and yu''^-' is the next-to-leading correction to the Aq self-energy. Note that Ui^2 will receive 
corrections from higher-loop diagrams. 

In [32], the authors computed (-R^)/? numerically. To match the notation there, we have 
Ui = Xi, and U2 = 3"'^'^ -|- |, where Xi = 2.298, Xs = 0.719 and Xa = —0.082 according to 
Monte Carlo method, leading to Ui = 2.298 and U2 = 1.818 in 



7 Further Discussions 

7.1 More Supersymmetric Gauges 

As discussed in [231 EH EH] , manifest off-shell SUSY is essential for solving the Schwinger- 
Dyson equations. The importance of having an off-shell description can be seen from the 



fact that our final solutions for the nonzero modes (1.6) respect the SUSY Ward identity 



(2.16), which only makes sense for an off-shell SUSY formulation, to leading order in large 

/3. 

A related fact is that having kept the auxiliary field / explicit, all the couplings in the 
A^ = 4 Wess-Zumino quantum mechanics are cubic. That is, the quartic interaction 0^ is 
represented by the cubic coupling /0^. It appears that having cubic couplings only is the 



key for the soft collinear approximation (3.1) to work. 

If we want to apply the soft collinear ansatz to Schwinger-Dyson equations for the BFSS 
matrix model and thus solve them analytically, it seems crucial to trade every quartic inter- 
action with cubic couplings by "integrating in" auxiliary fields. Recall that the BFSS action 
can be decomposed into a A/" = 4 vector multiplet V and three TV = 4 chiral superfields $" 



with superpotential W = — i-^g'yMe^^'^Tr ($'*[$'', $'^]), see (5.15). The quartic interactions 
among 0's are already traded with cubic interactions involving the auxiliary fields /. How- 
ever, there are still quartic interaction between the vector and chiral multiplets, e.g. from 
D^(j)"-D^^(f)"-. It would be desirable to trade these quartic couplings with cubic couplings as 
we did for the self-interaction in the chiral multiplet. 
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Manifest A/" = 8 supersymmetric formulation using the harmonic superspace [33] pro- 
vides a natural solution. In the JV = 8 language, the BFSS matrix model consists of one 
vector multiplet and one hypermultiplet. The action for the (gauged) hypermultiplet takes 
a particularly simple form, 



Shyper = " f dudC^-'^q+{D++ + V++)q+ . (7.1) 



g"*" is the (analytic) superfield representing the hypermultiplet and the + superscript labels 
a f/(l) charge which is not important here. D~^^ is a derivative acting on the w-space and 
dudC^~'^^ is the measure on the analytic superspace (the analog of chiral superspace in the 
Af = 4 case). ~ is a special conjugation whose details are again not important. At last, V~^~^ 
is the superfield for the vector multiplet. 

The key observation here is that, in the harmonic superspace formulation, the couplings 
between the hypermultiplet g"*" and the vector multiplet V^~^ are cubic. Note, however, that 
vector multiplet self-interactions involve quartic and higher couplings. 

The complication about the harmonic superspace formulation is that each superfield 
contains infinitely many component fields when expanded with respect to u. Therefore, to 
preserve manifest A/" = 8 SUSY, we necessarily have to deal with infinitely many coupled 
Schwinger-Dyson equations. An Af = 8 supersymmetric gauge fixing condition 

D++1/++ = (7.2) 

will kill all but a finite number of component fields in V^^. However, the hypermultiplet 
superfield g"*" still contains infinitely many component fields. 

It may also be possible to formulate Schwinger-Dyson equations for BFSS matrix theory 
with all A/" = 16 supersymmetries manifest, by working with a gauge fixed version of the 
Batalin-Vilkovisky action written in terms of pure spinor superfields [3l]. These are left to 
future investigation. 

7.2 Schwinger-Dyson Equations for Supersymmetric Field Theo- 
ries in Other Dimensions 

As another potential application of the technology we developed, one may employ the loop 
truncated Schwinger-Dyson equations to supersymmetric field theories in more than one 
dimension, and compute correlation functions and free energy along an entire RG flow tra- 
jectory. We illustrate this with the example of the two-dimensional A/ = (2, 2) Wess-Zumino 
model with $^ superpotential. 
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It is well known that the M = (2, 2) Landau-Ginzburg model, with the superspace action 

S= [ £xd'9 $$ - / d^xd'e ^^ - / d^xd'e ^^, (7.3) 

J J k + 2 J k + 2' 

flows to the k-th M = 2 minimal model in the infrared [5Sl ES]- The equivalence between 
the IR flxed point of the Landau-Ginzburg model and the Af = 2 minimal model has been 
checked in many ways: the central charge of the SCFT at the flxed point was derived in [36] , 
the elliptic genus of the two theories were shown to agree [37], the integrable RG flow of the 
2d A^ = 2 LG model was worked out in [3H], and the exact S'-matrices and the central charge 
were determined along the entire RG trajectory [SH]. However, apart from numerical results 
obtained via lattice simulation |10l E] , the understanding of this RG flow is based on either 
symmetry arguments (exact scaling dimension of $ in the IR limit from the R-charge [36] ) 
or very special integrable structure of the theory at the IR flxed point. The computation of 
correlation functions along the entire RG flow is not straightforward in such approaches. 

The loop truncated Schwinger-Dyson equations for the M = 2 LG model, on the other 
hand, appear to provide a useful approximation and expansion scheme along the entire 
flow. In particular, the one-loop Schwinger-Dyson equation already gives the correct scaling 
dimension of the flelds in the IR limit. For the moment, let us consider the zero temperature 
limit of the one-loop truncated Schwinger-Dyson equations for the self-energies of the matter 
multiplet. Unlike in the case of supersymmetric quantum mechanics, here the solution to the 
S-D equations is non-singular in the zero temperature limit. Clearly, SUSY Ward identity 



(C.24) is obeyed for the zero temperature solution, and we can express the fermion and 
auxiliary fleld self-energies in terms of that of the scalar fleld. The one-loop S-D equation 
then takes the simple form 

"*"' -fl^ ., \... ^ (T.4) 



p^ " J (2t,Y (q'- + o(q))((p-qf + a(p-q)) 
In the p — )■ limit, the equation reduces to 

a{p) ^g2 ff3_ __i (7.5) 



p'^ J (27r)2 cr{q)a{p — q) 

Hence we flnd the solution for the self-energy in the IR limit, 

^(^) ^ ^i/3r(|) 3'P'- ^^-^^ 

Indeed, we have rediscovered the | scaling dimension of $ in the IR. While it may seem that 
this simple computation of the IR scaling dimension is essentially equivalent to the statement 
that the cubic superpotential is marginal in the IR, the Schwinger-Dyson equation is not 
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specialized to the IR fixed point. One may solve (7.4) at all momenta p, which gives an 
approximation for the self-energy of $ along the entire RG flow. It appears that, by taking 
into account higher loop corrections (the two-loop correction is absent and the first correction 
comes in at three-loop), we will be able to obtain an expansion for exact correlators along 
the RG flow. 

In addition, one can calculate the free energy by solving the Schwinger- Dyson equations 
at finite temperature. In this case, we discretize the momentum p — )■ (p, ^). The S-D 
equations become an infinite set of coupled integrable equations. In the low temperature 
limit, scale invariance determines the free energy up to an overall constant, F = —ncVT'^/G. 
Here c is the central charge at the IR fixed point. It would be interesting to compute the 
value of c by solving the loop truncated S-D equations, and compare the result to the k = 1 
J\f = 2 minimal model value c = 1 (c = -^ for the fc-th minimal model). It would also 
be interesting to see how the two-point function for $ flows using the Schwinger-Dyson 
equations. We leave these for future development. 
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A Hawking Decay Rate of DO-Branes 

As argued in the introduction, a DO black hole decays by emission of DO-branes. The fact 



that these emitted DO-branes are unbounded by a DO black hole background (1.1) is due to 
the cancellation between the gravitational potential and the Ramond-Ramond potential in 
the DBI action for DO-branes. The decay rate of the black hole is governed by Hawking's 
formula [l2] 



n=l 



d^k 



V / -^^ -^^ (AV 



(27r)9 g^(w-M«) 



where k is the spatial momentum, u the energy, v = du/dk the speed (given by the dispersion 
relation at infinity), and n the DO-brane charge; a^hs is the absorption cross section and /i 
is the chemical potential for the DO-brane charge. The phase space and thermal factors 
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make the spectrum of the Hawking emission strongly peaked at a; ~ 1//3. The absorption 
cross section in the regime of interest can be obtained through a classical computation of the 
capture cross section, or a semi-classical (WKB) computation of the greybody factor and 
then summing over all partial waves. Here we present the latter. 

A.l Greybody Factor 



Consider a scalar field propagating in the DO black hole background (1.1). Because super- 
gravity can only be trusted for tq ^ N^'^lp, and this background is only a stable solution 
To ^ N^^^lp, we shall take N^^^lp <^ tq <^ N^^^lp]^ To incorporate the effects of the gravi- 
ton, the dilaton, and the Ramond-Ramond 1-form background altogether, it is convenient 



to lift ([LI]) to M-theory: 

fdx 
/ = „o 7 , Co = — , A 



2 
IIA 

?o - (1 + A)dxiodt + ^^"/^^ f-^dt"^ + A-^dr'^ + r'^d^l, 



(A.2) 



CqNPp . _ 15 



no' 



where xiq is periodic under xiq ~ xio + 27r_Rio, -Rio = QsU, and Ip = {27rgsY^^ls- Note that 
we have written all quantities in the M-theory scales Ip and Riq. To simplify the notation, 
let us set the Planck length Ip to 1, and rescale 

xio -> xioi?io, t -> t/Rio, cqN -^ N. (A.3) 



The Hawking temperature (1.2) is then 

1 7 fr^ 

^ 13 Airy N 
In this background, the Klein-Gordon equation for the n-th Kaluza-Klein mode is 

Ar^^dyAdr -2nuj(l + —] + -^-j = 0, (A.5) 



where L^ = i{i + 1) is the total angular momentum and 



R^^'^d-^] . (A.6) 



2n \ 2Nuj 



^^ When < rg < iV^/^, the solution (1.1) is unstable against the inhomogeneous IID black hole back- 
ground _5J. 
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Under a redefinition (p = Vr^Acp, the Klein-Gordon equation turns into a standard wave-in- 
an-effective potential problem 

-dj4) + [/0 = 0, (A.7) 

where 

U = A-^l\-^dl{r^A^I'') - 2nuA-^ (i + ^\+ ^H^. (A.8) 



ly ( J cf A 



The greybody factor is then just the transmission amplitude T of an incoming plane wave 
in this potential. Due to the complicated form of [/, we cannot solve the wave equation 



(A.7) exactly; instead, we will solve it for large and small r using the WKB method, and 
then match the solutions in an overlapping region. We show the existence of such a region 
presently. 



The potential \J (A.8) involves two natural scales, tq and R. When w <^ nrJ^/N, 

14 7 7 

and when u ^ nrl/2N, 



. Nu rl2Nu 7 ,, ^^, 

'^ =^ ar = T7i^ » '■»■ '*-i°' 

In either case, R ^ tq, so if we define the far region to be r ^ tq and the near region 
to be r ^ R, then we have an overlapping region ro ^ r ^ i?. The only part of the 

7 7 

spectrum where this fails is when w ~ |^. However, for a; ~ |^ ^ 1//3, the thermal 

7 

factor suppresses the emission of DO-branes exponentially, so the contribution from u ~ ^ 
is negligible. Similarly, the phase space factor suppresses the contribution from co < ^^, so 
we may neglect this region and assume nuR^ ^ 1. We will examine these two regions more 
carefully at the end of the section. 

In the far region r ^ rg, we can set A ~ 1. In terms of p = r /R, the wave equation 



(A.7) becomes 

-9^0 + ^0 = 0, (A.ll) 

where 

V = i^±^ - 2nujR^ fl + ^V (A.12) 



p2 \^ p7 



To check if WKB is valid, let us compute 



V.._2,11±^UM^. ,A.13) 



p3 J p8 

For L^ ~ 1, each of V and dpV is dominated by the second term since nuR"^ ^ 1, and so 
the WKB criterion |c)pl^p/|l^p is small for all p. Hence WKB is valid, and the far solution 
is 

0™ oc p-\-V)-'/'e^'^'P^^. (A.14) 

41 



In the near region tq < r <^ R, let us define p' = t/tq. Then the wave equation (A. 7) 
becomes 

-9j0 + Vr0 = O, (A. 15) 

where 

1 + 4(24 + L^)p'' + jSnuR^rDp'' - 4(12 + L^)p''' 

If L^ ~ 1, then the p'^ term dominates over 1 and the p'^ term because both R/ro and nuR"^ 
are large, and because the range of p' is bounded below by 1; the p'^ term also dominates 
over the p'^^ term because p' <^ R/tq. Thus the potential can be approximated by 

2nuR'^ p'^ ,, ^^, 

^- — n ,r,2 - A.17 

rl (1 - p'7)2 

The WKB criterion 

\d,W\^ ^^ rl 49(l + p'^)^_ ^^^^g^ 



1-^ 



\W\^ 2nujR'' p'9 
is peaked near p' ~ 1, the value thereof is small since both R/tq and nujR^ are large. Hence 
WKB is valid, and the near solution is 



'■■> 



0™ oc p'-4(-i^)-i/4e±^/'^p'v^. (A.19) 

Now we match the two solutions in the overlapping region ro ^ r ^ i?, which corresponds 



to p ^ 1 and p' ^ 1. But the two sides match trivially, since the two solutions (A. 14) and 



(A.19) become the same in this region. The WKB method can give a nontrivial reflection 



amplitude only if there exist turning points in the potential. Under our assumptions, the 



potential (A.8) is well below for the entire range of r, and therefore the transmission 
amplitude T is just 1. However, the potential (A.8) has turning points when L^ is large. 



From the approximate potential in the far region (A.12), this happens when 

12 + L^ > ( - j LLx, ^Lx = '^nuR\ (A.20) 



It is not hard to see that the exact potential (A.8) gives the same bound up to multiplication 
by an 0{1) constant and subleading order corrections. If L^ exceeds the threshold, then the 
potential barrier will provide an exponential suppression of the transmission amplitude. 



A. 2 Partial Wave Summation 

The absorption cross section is related to the greybody factor by ^3] 

oo 
CTabs = Xl^abs, (A.21) 
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where 



c^L = 2^7r^/^r(7/2)(^ + 7/2) (^ ) ^) J (A.22) 



is the absorption cross section for each partial wave. Let ^max be defined such that -^^ax = 
^max(^max + 1), then becausc T — I ioY £ < £max and T is exponentially suppressed for 



> ^max, we can approximate the sum (A. 21) by a hard cutoff 

^max 2^7r^/2r(7/2) ""^ 

'^abs — 

where 



aabs-2^a,,,_ g^g, ^g , (A.2dj 



Cax ^ /^ ^ v/2;^i?. (A.24) 

Note that even though we have kept the prefactor, corrections to the hard cutoff approxi- 
mation will multiply this prefactor by an 0{1) constant, so one should not take its value too 
seriously. 



A. 3 Decay Rate 

We are now ready to compute the decay rate. The dispersion relation at infinity is w = k'^/2n, 



so V = k/n = ^y2uJ/n and d^k = ((2n)^/^7r^/^/r(9/2)) co'^^'^du. The chemical potential is 



12 



7 
f,dt = C^\r=ro = -^dt. (A.25) 



Putting these together, the decay rate is 

r^^E^'^""'^°^'''/ ^^4^- (A.26) 

n=l " 

Due to the exponential suppression from the thermal factor, the dominant contribution to 



r comes from oj ^ nrl/2N, where (A. 6) can be approximated by 

,14 



nrn 



R' ~ ^^-^ (A.27) 



^^ It is the difference between the gauge potential at the horizon and that at the "far-horizon" infinity, 
which is by our gauge choice. Alternatively, without going to far-horizon, it can be fixed by the first law 
where one uses the regularized ADM mass for the black hole. 



43 



Then 

o4/7 ^16 °° . r°^ , ,20/7 



o4/7 ^16 ^ /"oo 

8!7r A^8/7 ^ Jo 



n=l 



^ 2^/^r(27/7)C(27/7) r^^ ,,./,^ ^^.28) 

Slvr iV8/7/327/7 ^ ^ 

^ 2^/^r(27/7)C(27/7) / 7 ^^^359/14 ^^^^^^^^^_ 
8!7r yATx) Ar43/i4 

where the sum over n is replaced by the n = 1 term because /3rl/N ~ yr^jN ^1. To 



conform with the units in the rest of the paper, one can undo the rescahngs (A. 3 ) and restore 
Ip. 

Now let us examine the regions u ~ nrl/2N and u < l/nR"^ where the computation of 



T in Appendix A.l fails. For u ~ nrl/2N, we lose the overlapping region because the two 



terms in the brackets in (A. 6) may cancel to give a small R. However, looking at the range 
of L^ (A. 20) for which T becomes exponentially suppressed, we see that a smaller value of 
R only lowers the hard cutoff and makes the contribution to F smaller. Since we already get 
an exponential suppression from the thermal factor, we can safely neglect the contribution 
from this region. For u < 1/ni?^, the cutoff is at most 0{1), so the suppression from the 
phase space factor again allows us to neglect its contribution. 

B ID A/" = 2 SUSY 

In this appendix we will discuss the structure of ID A/" = 2 SUSY, the vector and the matter 
multiplet. It is previously discussed in |23]. We will work in the Euclidean signature. 



B.l Superspace 

The (Euclidean) ID A/" = 2 superspace consists of one bosonic coordinate r and two real 
Grassmann coordinates 6a in the vector representation of R-symmetry group S'0(2)/j. We 
choose the S0{2)fi Dirac matrices to be real and symmetric, 7^ = cxi and 7^ = a^, where a* 
are the standard Pauli matrices. 

There are two invariant symbols, 6a/3 and e^/j, which represent the inner and exterior 
products for 2D vectors. We normalize them by 611 = 622 = 1,^12 = ^21 = 0, and ei2 = 
—£21 = l,eii = £22 = 0. If we were to distinguish upper with lower indices, there is a 
potential confusion whether one uses Sap or e^/j to raise and lower the index. To avoid this 
ambiguity, we will write every index downstairs. 
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Let us define 

so that 6'^ is realj^ We define our fermionic integral measure such that 

£9 e^ = 1. (B.2) 

Given a 2 x 2 matrix B^is, we can expand it as 



1r ^ +^ 

Some useful identities are 



-Ba/3 = -^Sppdap + ^ll^Bpalali + -^d-puB p„ea(5. (B.3) 



7a/37pa = ^apSpo- + ^aa^lSp " (^a/Jf^pa, 



The supercharges and super-derivatives are 



They obey the algebra 



(B.4) 



Ct7a 



{Qa,Q4 = -25«A, (B.7) 

{Q^,Dis} = 0, (B.8) 

{D,,D;3} = 25„^a.. (B.9) 



B.2 Vector Multiplet 

B.2.1 Connections 

The vector multiplet is constructed as follows. We first introduce a real connection on 
superspace 



Va = D^ + T^, (B.IO) 

\/^ = dr + iT^. (B.ll) 



^■^By real we mean the it is real when we Wick rotate to the Lorentzian signature. This will always be 
assumed implicitly in the following. 
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Note that F^ is an anticoinmuting real superfield while F,- is a commuting real superfield. 

For reason that will become clear in a minute, F^ is determined by F^ by a constraint 
equation, so we only have to expand the real connection F^ as 

F„ = xa + iAo9a + X'l'o.pOp + dea^Op + 2eap\p9\ (B. 12) 

The super gauge transformations are 

5fj:^ = iVak = iDc,k + i[T^,K], (B.13) 

5tJ:r = Vr^ = ^r^ + i[Vr,^], (B.14) 

with A a real superfield. If we expand 

K = Vl + iujJc. + We'^, (B.15) 

then the super gauge transformation in the component form is 

6AXa = ^(x + i[Xa,^], (B.16) 

5AAo = r2 + i[Ao,fi] + ^{Xa,Wa}, (B.17) 

5aX^ = i[X\n] - \^'^p{Xa.ujp}, (B.18) 

S^d= ~W + i[d,n]- -e^piXa.ojp}, (B.19) 

(^aA^ = -2^a - 2ea/3[X/3, VT] +i[A„,fi] + -[Ao,Wa] + -i'^p[X\ujp] - -e„;3 [rf, w^] . (B.20) 

B.2.2 Supersymmetry Transformations 

The SUSY transformations are generated by Qa in the following way 

8e^a = ^/jQ/jFq 
=> 6eXa = i£aAo + s p^^^X' + Speapd, (B.21) 

4^0 = e^Aa - -SaXa, (B.22) 

1 . 

6eX' = -i7a/3^aA/3 + -7a/3^«X/3, (B.23) 

5ed = -ieapEaXp - -ea[s£aXi3, (B.24) 



1 



(^sAq, = -e^Ao + -epP-i'^p - -epdeo^p- (B.25) 
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B.2.3 Field Strengths 

The field strengths Fajs and J-^ are 



{V«, V/3} = 25,^ V, + J-„/3, (B.26) 

[v., V„] = -iJ'a. (B.27) 



Explicitly, they are given by 



J-„/3 = Dj:p + DpT^ + {r„, r^ - 2i5^pT^, (B.28) 

Fa = tdrTa + D^Tr - [Tr, r„]. (B.29) 

Under super gauge transformations, the field strengths transform covariantly: 

5A^a/3 = ^[^a/3,A], (B.30) 

5aJ'« = «[J'„,A]. (B.31) 

Note that the lowest component {6 = 0) of the field strength is neutral under the super 
gauge transformations which depend on 6. 

Next we impose the constraint 6apJ^ap = 0, from which we can express F^ in terms of 
others 

r. = -^(Z},r, + r,r„). (B.32) 



Note the super gauge transformations (B.13) and (B.14) are compatible with the constraint 



(B.32). Now J^ai3 is a symmetric traceless matrix, so it can be viewed as a vector in S0{2)ji 



by defining J-"* as 



J"' ^ h'apJ'a^ = llap (^(^r^) + r(,r^)) (B.33) 



B.2.4 Action 

The ID A/" = 2 super Yang-Mills action is given by 



Sym =^ jdrd'eTi I - ^ (^e^^ V^rV^r] (B.34) 

In the Wess-Zumino gauge, the component fields Xa and d in Fq are gauged away by ap- 
propriately choosing the super gauge parameter Ua and W. While the part of the super 
gauge symmetry is fixed by the Wess-Zumino gauge, we still have the conventional bosonic 
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gauge transformation generated by Q. In the Wess-Zumino gauge, the expansion of the real 
connection superfield reduces to 



The field strength J-"* in the component form is then 

The Yang-Mills action in the WZ gauge is 

Sym 



9ym 



drTi { ^KDoK - ^\alL^[X\ \p] 



+ ^DoX'DoX' - ^[X\X^][X\X^] 



(B.35) 
(B.36) 

(B.37) 
fB.38) 



where Dq = d^ + i[AQ, ■ ]. Note again that (B.38) is only invariant under the conventional 
(bosonic) gauge transformation, but not under the super gauge transformation. 

Next we want to write down the action without fixing a gauge. Define the "covariant" 
component fields to be 

1 



X^lcov =J^^\e=0 — 7^1al3 {D{a^i3) + r(Q,r^)) 



61=0 



--^' + ^'^apXaXl^, 



A, 



O\cov 



=Ao 



e=o — —i-^ {DaTa + TaTc 



"^ nXaXo 



Aq 



=-F, 



Q 61=0 



idrTa + DaTr - [Tr,Ta\\ 



e=o 



2 1 % ■ ■ Z % 



(B.39) 



(B.40) 



(B.41^ 



To obtain the vector multiplet action without fixing a gauge, one has to find a extension 



of (B.38) such that it is invariant under the full super gauge transformation. Suppose we 
replace every X*, Aq, and Xa in (B.38) by X*|coi,, ^o|co»;, and Xa\cov defined above, then 



the action is still invariant under the conventional gauge transformation because the gauge 
indices are properly contracted and the derivatives always appear in the form of gauge 
covariant derivatives Dq. In addition, since the covariant fields are the lowest components of 
J-"*, J-'q, and r,-, they are neutral under gauge transformations that depend on 6. As a result, 
the new action is invariant under the full super gauge transformation as desired. Thus the 



ID A/" = 2 super Yang-Mills action is given by (B.38) with every component field replaced 



by the covariant component fields defined in (B.39 B.41) 
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B.3 Matter Multiplet 

Consider a real scalar superfield $ in the adjoint representation of the gauge group: 

^ = (j) + i%l)J^ + ife'^. (B.42) 

where the i in front of / is such that the kinetic term for / has the right sign. 
The SUSY transformations are 

4$ = e^Q^^ (B.43) 

^ 5,0 = -i£«Va, (B.44) 

^ei'a = -iepfe^p + 16^,(1), (B.45) 

4/ = £aea/3V'/3- (B.46) 

The super gauge transformation is 

(5A$ = i[$,A] (B.47) 

Written in terms of the component fields, we have 



5A0 = i[0, f^], 


(B.48) 


^Ki^a =i[V'a,fi] +i[0,Wa], 


(B.49) 


(^a/ = ea/3{V'a, W/?} + ^[/, ^] + ^[0, W^' • 


(B.50) 


The action we will consider is 




./ V 4 3 / 


(B.51) 



where a is in some representation of the flavor symmetry G and e'*^'^ is a totally antisymmetric 
G-invariant tensor. If we write the BFSS action in the M = 2 language, then it consists of 
one vector multiplet and the matter multiplet in the 7 of the the flavor symmetry G2. 

B.4 Schwinger-Dyson Equations for BFSS in the J\f = 2 Language 

Let us write down the one-loop truncated Schwinger-Dyson equations for BFSS with manifest 



A/" = 2 supersymmetry. The full action is ^vec + 5'mattcr + •S'ghost given in (5.1 ) and (5.6). Let 
us denote the propagators for the vector, matter, and ghost multiplets by 



vector : o"^, a^ 



matter : A^, g^, el, (B.52) 

ghost : Sn, tr, Un, 
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and write them in terms of the self-energies as 

1 



a' 



A 2 9 

" ^ pTrnp I J^' 5''' ^ ~ 27rr | yM' ^n ^ ]^ _^ ;=;M ' (B.53) 

1 i 1 



There is also a propagator for the zero mode of the gauge field, (Aqo^oo) = Po = Va^j where 
Aqo is normalized by Aq = P~^Aoo. 

The one-loop truncated Schwinger-Dyson equations are 
n-=2^ 1 3^ 1 14^ 1 



14 ^ 1 1 ^ 1 

2 8 1 

y _ 3 27rr ^ 1 4 27rr >^ 1 14 27rr ^ 



/3 /3 t'(^)2 + nv; /3 /3 ^(2|£)2 + _^sr /3 /3 tr(^)2 + n^ 

5 27rr 8 27rr 1 14^ 1 

14 ^ 1 1 ^ 1 



-y- 



1 



(B.54) 
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n ^ Z^ (2^)2 + n^ /3 4- (^)2 + £j:V ^ /3 a. ((2^)2 + n^)(l + HM 

12 v^ 1 4,2?™ 



2 8 1 

+ ^~^i + (^)2nM' 



2 ^ 1 4^ 



1 



^M _ 6 -^ 1 + - V 



(B.55) 

c_1.27m2V- 1 2 1 

«l /P ^ Z^ 



/5 /3 ^ic-fr+^rr-^+^^^r) ^f^^+a.f^^" 



■^G 



1 ^ 1 ^>r- 



1 1 (B.56) 



^2/3^fl 



1 



+ (oi)^n^)((?^M)3 + sr_„) 2/3/. ^ + Ef 



/O Z^ 



^ 2TTn 

1 



/i_2^ 1 _5v ^ '^V "^ 

2 " ^^ (W+nT " 2^ V (W+^P? ^ ^ n (¥)^+nF 

4 ^ 1 4 ^ 1 14 ^ 1 

R 2-^ {2-Kn _|_ /3 TT^'lS /? -^ |'27rr , / /3 ^2vy^2 P, 2-^ ("inn , /3 YIM\2 (-D-57) 

7^1 1 ^ 1 1^1 

+ ^2^(2^ + SM)2+^Z.(2^ + _|_nG)2-^Z.(2^ + SG)2- 
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B.5 Free energy of BFSS in the J\f = 2 gauge-fixing condition 

The free energy can be approximated using the mean- field method [23] . It consists of four 
parts, 



/3Fo = - -log^-^logaf + ^loga^- -^logA^ + r^log^^ 

I r I r 

7 V 



I Ij^O r I 



r,s l,r l,m l,m, 

14i v^ « 9 27rr 7 v-^ A? 5i v-^ Svrr a^ 2 v-^ (^f 



' I ' ^ 1-^ r ' ^ ' l+r+s=0 ^ 

^ l+r+s=Q ^ l+r+s=0 ^ I ^ ^ 



4 



1 V^ f 2 /2vr/ 2 \ 

+ - 2_^ \ai tr-ts - Ultras + [—) SiUas I 

^ l+r+s=0 ^ P / 



(B.58) 



(B.59) 



(B.60) 



|E| + f E A?A^4-f E A?*9- (B.61) 



2/3 ^/i ' /3^' P ' /i' 

Here /3Fo is the free energy of the trial action, PF2 is that from the difference between the 
quadratic part of the BFSS action and the trial action, and PF3 and /3F4 from the cubic and 
quartic couplings. 
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C ID A/" = 4 SUSY 

The ID A/" = 4 superspace is the same as the 4D A/" = 1 superspace, so we will mostly adopt 
the 4D notations here. We will work in the Euclidean signature. 

C.l Superspace and Convention 

Let us start with the 4D M = 1 superspace. We will denote (2, 1) and (1, 2) of 4D rotation 
group 5*0(4) = SU{2)i x SU{2)r by the lower a and a, respectively. The 5*0(4) invariant 
symbols e^/j, e^6 and their inverses are normalized by 



6^2 = e'^ 



-1, ei2 = ei2 = -l- (C.l) 



We also have two other 5*0(4) invariant symbols cr^^ and a^"" with /x = 0, 1, 2, 3 defined as 

'"' ^ '''f; (C.2) 

Under dimension reduction, the time direction is singled out and the 5*0(4) rotation group 
is broken to the 5*0(3)/? = 5*^/(2)^ i?-symmetry of the ID quantum mechanics. In addition, 
the 4D A/" = 1 SUSY theory has an t/(l)ij symmetry, so the full /^-symmetry of the ID 
A/" = 4 quantum mechanics is f/(l)i? x 5*f/(2)^. Note that the SU{2)fi is embedded in the 
4D rotation group SU{2)e x SU{2)r in such a way that the zeroth component of a^ and a^, 
i.e. i(5^a, is its invariant symbol. 

The supercharges Qa^Qa and the supercovariant derivatives V^^Va are defined as 

f^ (C.3) 

The only nontrivial anticommutators are 

{Qa,Qa} = 25^^dr, {V^,V^} = -25aadr. (C.4) 

As opposed to the M = 2 case, we will distinguish the upper index with the lower one. 
The indices are raised and lowered by the invariant symbols e^/j, e"-^ and their inverses, 

t/'" = e"%, V^a = ea/3^^ (C.5) 
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and similarly for the dotted indices. We introduce the following notation for the spinor 
bilinear, 

X'^ = X^'tpa = -e°''^Xaipfs = V'X, 



In particular, 



V^x = 


- ^aT = e^'^^aX^ - 


= x^- 


9a9i3 = 


-66eai3, d°'6 = 


-Uee-f', 


Oadp = 


-\oKp. 0^0^- 


Zi 



(C.6) 



(C.7) 



C.2 Chiral Multiple! 

The chiral superfield $(r, 6*, ff) satisfies 

P,$(r,e,^) = 0. (C.8) 

We can solve this constraint by first introducing 

2/ = r + (5,^rr. (C.9) 

By using T>cfioL = and "D^y = we find that any superfield $(?/, 6) that is a function of y 
and 6 is chiral. 

We can expand $(y, 6) in terms of component fields, 

^^{y, 6) = riy) + v^<(z/)^" + ^ny)o^ (CIO) 

with a being the index of some flavor symmetry group G. Note that we are working in the 
Euclidean signature so there is an i for the auxiliary field f{y) to make its quadratic term 
in the action positive. 

The SUSY transformations are 

[Q„, r] = ^v^c {Qa, r^} = V2e^pr, [q„ n = o, 

[Q^,0i = 0, {Qa,rp} = -iV25p^r. [Qa.n = -v^w^ 
[Q^, r] = ^v^c, {Qa, v^j} = - v^e .^r , [Qa, n = o. 
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The Lagrangian with a general superpotential Vr(<l') is given b}^ 

c= j (feSe <!$ + ( / Se iy($) + h.c. ) . (c.12) 



^^Note that the complex conjugate is defined by first Wick rotating to Lorentzian signature, taking the 
complex conjugate, and then Wick rotating back to the Euclidean signature. 
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C.3 Vector Multiplet 

Consider a real superfield V{t, 6, 9) in the adjoint representation of some gauge group with 
the following component field expansion 

v(r, e, 9) = r{t) + ex + Oxir) + eeM{T) + emir) 

1 - (C.13) 

+ r^^ (iS^hM^) + <6^'(^)) + dddK^) + Oee\{T) + -9999D{r). 

The gauge transformation is 

where A and A are chiral and anti-chiral superfield, respectively. 
The field-strength superfield is defined by 

W^ = - I'D^V'^e^^V^e-^^ (C. 15) 

8 

The kinetic action for the vector multiplet is then given by 

C = Ti ( f d^9 -^WWa + h.cA . (C.16) 

C.4 A/" = 4 SUSY Ward Identities 

The A/" = 4 SUSY Ward identities for the chiral multiplet can be similarly derived as in the 
A/" = 2 case. We define the self-energies o"(p), ?7(p), and /i(p) for the boson 0", the auxiliary 
field /", and the fermion ■j/'^, ■?/'"° to be 

p^ + (T(p) 
rab 
inp)f\-p)) - ^^^y (C.17) 



(C(p)^J(-p)) = -I 



5"^, 



al3 



p + hip) ' 

respectively. Note that o-(— p) = o'(p), rj{~p) = rj{p), and h{—p) = —h{p). 

The SUSY Ward identity is an exact relation between the exact propagators for fields in 
the same SUSY multiplet. Suppose the SUSY is unbroken in the model, we have 

O=({g.,0'^(r)V5J(r')})=^v^(C(r)V^^(r'))-^v^5„,(r(r)0^(r')). (CIS) 
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Going to the momentum space, this imphes 

(vsj(p)C(-p)) = -^ps^p{r{p)<P\-p)) (C.19) 

that is 

a{p)=ph{p). (C.20) 

This is an exact relation between the self-energies for the boson 0" and the fermion ip'^. 

Similarly, we can consider 



In the momentum space it is 
Using e^P5a'^5pp = —ea/3, we have 



v{p) 



Hp) 



P 



In summary, we have obtained the exact relation 

hip) = =Pvip), 

P 



(C.2i: 



(C.22) 



(C.23) 



(C.24) 



which takes exactly the same form as in the A/" = 2 case (2.16) 



C.5 Vector Multiple! Zero Mode Action under J\f = 4 Gauge- 
Fixing 

In this subsection we will truncate the vector multiplet action to the zero mode sector in the 



A/" = 4 gauge-fixing condition (5.16). We only need to focus on the component fields Aq,X^ 



and R since D = Aq — ^R is set to zero in the zero mode truncation. We will write Aq and 
X* collectively as A^ with /x = 0, 1, 2, 3. 

The field-strength superfield is 



Wa = --V^V^e^^Vo.e-'^^ . 



(C.25) 



To get the terms with only zero modes, we can replace every V in (C.25) by d/dO and 
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similarly for V. Using 
exp (2R + 2e''e^a^^■A^ 



e2R+ /"'rfse2^^2r^VM p2(i-s)« 
Jo "^ 



+ f ds f rfMe2"^2r^^a^A^e2(^-")^2^^^V:;^A^e2(i-^)^ (C.26) 

•/ U •-' u 



e^^ + 2 / rfs e^^^ r^^a^ A„e^(i-^)« 







a/3--^"- 






we then expand 



Now 



• 1 rl 



(C.27) 



^ Jo io 

We can simplify the above equation by using the identity 

^'%<e = ^a.ay^'^ = -e^,5^^ - 2te^,{Sn,^ (C.28) 

where 



+ f rfse2^%e2(i-^)^ f due-^'^'A^^e-^^'-'^A (C.30) 

Jo Jo 

where we have used ea/se'^^ = —^^'^f^ + 5^5^ and {S^)^ = 0. Note that in the case of 
R = 0, only the term with S'^" survives and is contracted with [A^, Ay], as expected in the 
usual Wess-Zumino gauge. 
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The zero mode terms in the Lagrangian can then be written as 
1 



£ 



Ag: 



YM 



e'^^WoWfi 



+ h.c. 



where 



We have used 



9ym 



-2 / ds I rfMe2(i-")«A^e-2(^-")«A'^e-2(i-^)« 



'0 JO 

Jo Jo 



9ym 



J±^ / ds^±2sR, ±2[l-s)R_ 



with e^^P" normahzed as e" '^ = 1. 



0123 



(C.31) 



(C.32) 
(C.33) 



D Computation of Xn 

In the following we will frequently encounter divergent sums over integers or half-integers. 
Since these divergences come from expanding propagators in series of say ^r?, the sums 
should be understood as regularized sums cut off at n ^ a/Z^/ctq- 



Using (|3.5|), (|3.7) and (3.8), the equation for x„ is 
2 1 



[l + Xr. 



where 
A, 



Trn \ •'— ^ 



liao'^+Sn{l+Xn)-^ 



n — k 1 1 + X 



+ A 



/5c7o ^ + Sn{l + Xn) 



+ Bn 



(D.i: 



T™ \^^^ k SkSn-k 1 + Xfc 



+ Xn-k Sr~^ 1 



27r|r2|?7o 



13 (3ao 



[1 + o{p 



-6/5 



^,4;^ kSkSn-k 0- + Xk){l + Xn-k) 
k^O,n 



(D.2) 



Since A^ contains r/o, to solve for a;„ we also need to know tjq to the relevant order, 

^0 = -5^ + 0{P"% (D.3) 
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The Schwinger-Dyson equation for Xn can thus be written as 
2 ,27rn,(l + a;„)2-l 



/3^o' /3 



A„-(l + x„)25„ + 0(r''/'), 



(D.4) 



or 



/3 ^22 



'27rn' 13 



n 



. ,_ rCShSn—k \ J- 



k=/=0,n 



i + Xfi—k 



[1 +Xr. 



kSkSn-k \ l + Xk {I + Xk){l + Xn-k) 



E 



1 1 + Xn^ 



L^ h h 



hrhn-r '27c\n\r]o V 2 



/3 /3ao 



(D.5) 



Note that the correction in replacing s„ by s„ is of order (9(/3^^/^) because A„— (l+x„)^-Bn ~ 
/3-1. Write 

(D.6) 



Sfc = sign(A;) 



— + 4^ /.. = sign(r)^— + /.«, 



where s^. ,hr are of order (3 "^1^ . 

Sn can be determined by expanding the Schwinger-Dyson equation for Sn as written 
earher. Let us also write for now the Schwinger-Dyson equation for hr to the relevant order, 

2 1 1 



hr 



(3(Toh, 



ly^ 1 



kSkil + Xk)hr-k 



+ o{r'). 



(D.7) 



We have 



(1) _ (3cro ^^ sign (£) sign (/c - £) 



Air 



i^O,k 



+ o{p-'] 



h^}^ 



P(^o y^ sign(f)sign(r - i) ^ ^/^-l^ 



(D.8) 



47r 



1^0 



The equation for Xn can be written as 
(l + x„)'-l 
, /3 ^22 



•2?™^ (3 
_l3ao 

,(3cro 



/3ao >^-^ sign(A;)sign(n - A;) A + x^-k 



k^Q,n 



k 



^L -\- Xjij 



y sign(/c)sign(n — k) i — 



1+Xk {1 + Xk){l + Xn-k) 



+ Xk 



+2( ) 2 \^ sign(n — k)sl, — N^ sign(n — r)hl^' 



\k=/=0,n 



27r\n\r]o V 2 



/3 /3ao 



C(/3-3/5). 

(D.9) 
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Note that we have set Xn to zero in the third hne since it is already of order /?". Using 
2_] sign(n — A;)s^ — N^ sign(n — r)hl^' 



_ /3oo 
~ An 

we have 



^ \p\ 



e^O ' ' lky^O,n 



y sign(n — A;)sign(A; — i) — /^ sign 



27rV [(1 + x^r - 1] = ^n J2 



k^O,n 



k 



{n — r)sign(r - 


-i) 


+ c?(r'), 






(D.IO) 


1 + Xn^k 


(l+:r„)2 ^ 



l + ajfc (1 +a;fe)(l + x„_fc). 



+ 



+ 



,/3o-o,5/3 



TT 



y- 



y sign(n — A;)sign(A; — i) — 2_. sign! 

.k^O,n,t. r 



n — r Sign r — 



/3Vg //3ao 
27r|n| 



0(/3- 



-3/5N 



(D.ll) 
where we have used r]Q = 1//3^(To + 0{(3^^^). Note that on the RHS, even though the /3^(To 
term is naively of order (3^^^, but the two pieces in the bracket cancel to leading order, and 
so the term is actually of order /3°. Without further cancellation, the whole RHS is of order 
1. On the other hand, the LHS is of order (3~^^^. Therefore we have an inconsistency unless 
the RHS vanishes to leading order. 

Write 

a;„, = f^)i^ + 0(/3-6/^), (D.12) 



2 vr 



where a„ — )■ as n — )■ oo. We assume that a„ is of order 1. The equation for x„ becomes 

Esign(A;)sign(n — A;) , n 3 ,^,^ o/r, 

I i^n-k - «") + y + ^(r ^ )• (D.13) 

Here we used 



2^ \P 



e^o 



\i\ 



y sign(n — A;)sign(A; — £) — N^ sign(n — r)sig: 



.k^O,n/ 



1 1 

Y] Tji (sign(n)sign(£) + Sne) = ^• 
^-^ \t\ \n\ 



(D.14) 



£7^0 



ttn now obeys the equation 



sign(n)/JA^(|n| — 1) 



1 
2n 



sign(nj 



n^ 



2E 



sign(fc)afc 

\n — k\ 



(D.15) 
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where 



A solution is given by 



that is, 



m ^ 

HN{m) ^Y^l- 



k=l 



2\n\ 



xLr, 



,(3(To. 



+ o{p- 



-6/5^ 



(D.16) 



(D.17) 



(D.18) 



' 2 ' 27r|n| 
This confirms our intuition that SUSY Ward identity should hold for the nonzero modes to 



leading order in the low temperature expansion (2.17) (x„ essentially measures the violation 



of the the SUSY Ward identity). Together with (3.7), (D.6), and (D.8), we have solved a, 



rjn, and hr to the next-to-leading order in large {3 expansion in terms of o"o. 



E Computation of s« , ft« , sj? and h 



(2) 

r 



In this appendix we will compute s„ and hr to the second order, (9(/3~^). These terms are 
relevant for solving the Schwinger-Dyson equation for ctq. Let us collect the equations for Sn 
and hr to the relevant order. 



2 ^2 1 / 47m 1 



'/3ao si 
. 2 .o 1 



hr 



l3ao s3 
2 

2 1 



1- 



(3 Sn 

Ann 1 

/3 Sn 



^n^n j_ ^ ^ 



^n-Dn 



+ 



(3aosl 
2 1 



1 1 ^ 1 



kSkil + Xk)hr^k 



iSaosl 



+ Bnil + Xn) 
l+Xn 

(An + Bn) + 0{(3-^"') 



o{p 



-8/5N 



47rr 1 ^^ 



/3o"o Y^ sign(A;)sign(r — k) 



E 



/3o-o hr P'^a^hl n ^ kskK-k 27r ^_^^ 



A; 



a:fc + 0{(3 



-8/5^ 



(E.i: 



^n = — I ^ y^ 1 y^ — \- y^ i—r 



+ 



27r ^ fcra 



I Ly I Ly-t ■)-■ 



^;;:-sign(A;)sign(n - A;)(x„_,. - Xk) + _J, H ^)" + 0{p-^"'), 



^2j^2 



-B« 



Define 



E 



" k^.n ***"*—* 



5.0 j^ ^'i^"(*^)y("-'0 (,^,^,^_^,)^o(^-»/.) 



27r 



k=/=0,n 



k 



(3cro r^ _ ST sign (f) sign (A; 



(E.2) 



e^o,k 
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We can then write 



(E.3) 

Zn 



M ^^A ;,(2) 



with Sn and hn oi (3 ^ order. An and i?„ can be expanded as 



-'Ic - 

■n '^n 
TT 


27r2 ^ 


+ - > 

ky^O,n 


n — k . 
sigr 
kn 


> _ 9^r 

'n ^n 
TT 


27r2 4^ 



sign(n - A;)Cfc + sign(A;)C„_fc , fi- v^. -.xfc- / ^n^ 

; r— > —1 Sign n Cfe 

/7 \ • / 7 \/ \ 2sign(n) ^ .^, o, 
sign(A;)sign(n - k)[Xn-k - Xk) H ^^fi- + els' ), 

sign(n - A:)Cfc + sign(A:)Cn-fc 
k 



_g;^ sign(^)sign(n-^) ^^^^^^_^^^^^^3). 

fc/0,n 

The equations for Sn and /i^ are equivalent to 

47r72 

((7s 0^ = 1 - -— + g^AnBnsl + ^'s„(A„ + 5„) + O^g'^) 

P^n 

Aim , 5-^ „2 , '^9'^SnCn 9'^Sn >r^ sign(n - k)Ck + sign(A;)C, 



(E.4) 



n—k 



, ^^s„ v^ . ,xfc . / k^r^ 9*Sn sr^ sign(A;)sign(n - /c) 

+ ^— > (-1) sign(n - -)Ck > [Xk + a;„_fc) 

fi'^'sn v^ ^~^ . /7\ . / 7\/ \ 2sign(n) 7 ox 

H > , —. sign(/c)sign(n - k)[Xn-k - Xk) H t^^9 Sn + 0[g ) 

TT ^-^ kn T{^n^ 

ky^O,n 

^ ^ _ 47rg\n\ ^ 2ghign{n)Cn ^ ^g'^Cl 



(E.5) 



/3 TT 7r2 

9^ . / X v^ sign(n - k)Ck + sign(A;)C„_fc 
^sign(n) ^ 

k^Q,n 



, 9^ \^ / lAfc • / ^\n 2^^ • ^ ^ Y^ sign(A;)sign(n - A;) 
+ ~^ir~\ 2^ -^ '^iS'^ "- ~ o '^^ sign n > Xk 

' ' kj^0,2n ky^O,n 

+ ^g' + 0{g'), 
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2 T 27rr g^hr ^^ 1 9^ ■ f \^^ sign(A;)sign(r - k) g 



r . / X V^ sign(A;)sign(r - A;) g 

- — sign(r) 2^ Xk + 0[g ). 

(E.6) 






Sn and /if. can then be solved to be 

(2) fi'TT^ , ■ , ^ 9^ r<2 9^ ST sign(n - k)Ck + sign(A;)C„_fc 
^4 ^ = - ^ + sign(n)— Q - — )^ ^ 

, / Y^ ^ Nfc • / ^\n 9^ v^ sign(A;)sign(n - A;)xfc / g 

(2) 9T^r , . . ^ / ^2 / v;^ sign(r - A;)Cfc + sign(A;)C^_fc 

_ 9^Y^ sign(fc)sign(r - k)xk ^. g. 

^'^ MO ^ 

(E.7) 

It follows that 

/ v-/ r.kr. \^ sign(n - l)sign(n) 27r^ >^ ^£ 

- ^ Z.(-i) ^1 Z. ^ + -J- Z.(-i) Is' 

?' ^. .M- ..N ^ sign(f-A;)C, + sign(A;)C| 



, 5'^ V^ • / N V^ sign(A;)sign(n - k)xk g^ yr^ . / x v^ sign(A;)sign(r - A;) 

Uj^O kj^O,n r k^O 

- ^/ + Oig'^) 

= _i _ !1Z _ M 
2/3 6' 

(E.8) 
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where we have used 






1^0 



siga{k)Ci 
k 



Mo,f 






(E.9) 



F Computation of the Free Energy 



We will drop the overall N'^Nf factor in the following. Let us compute PFq first, the free 
energy for the trial action 5*0, 



/3Fo = ^ln^ + ;^ln(l + r7o) 



2 47r2 2 

oo 



n=l 



27rn 



P 



~r ^n\ -'- ~r '^r; 



+ Ei° 



n=l 



27rn 



/3 



~r ^nv-L ~r *^n 



,-1 



^ v^ , / 27rr 

2 2]ln( ^^ + /i, 

r>0 



/s 



co„=t + s ( 5: .s - 5: /.[;,' ) + 9^ ( E " - z -- ) + 4('3-«)' + o(fi-"% 



vJl^^O 



/3 



»?i=l r>0 



192 



(F.l) 

We leave the second and third term in the last line undone because they will be cancelled 
by the terms in —\{SI)q. 

To compute — 1(5*1)0, it will be convenient to recall our solutions for the propagators. 



A, 



/3 



2nn^-f + Sn{l + Xr. 



/3 



-^g 






27r|n 



27r|n| 



+ 0(/3- 



A(o) 

n 

A(2) 



/3 



27r\n 
27r\n 



■9 



Ad) 



/3 / ^V ^ 



27r n \ 27r 



f (2) , ^V 2^l^l , 9^ n^ , A 



(F.2) 
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27m 



/3 ^ + s„(l+x„) 



-1 



2n\n\ 



/3 



-9 



ff' ry , _(2)vT ^ , ^2^ . 27r|n| 



.(0) 



27i\n 



C\n\ + gsl){l - Xn + Xj + g- 



p 



-1 



+ 0(/3- 



(F.3) 



/3 



-^, e 



(1) 



5' TT^VA - a^n , 



/3 ^ V 27r 



27r|n| / (2) 5-^ 27r|n| , / 



y- ^ i^^^ini; • 



27rr 



+ /ir 



g signer j 



(2) 



1 + ^C\r\ + (7/lL + (? 



2n 



27r|r| 



-1 



^^°^+^^^)+^^^^ + W-^), 



gf'> =gs\gn{r-), g^^^ = - 



sign(r)— q^l, gf^ = g?>ign{r) {-gh^^^ + ^C'h - ^~^ ) • 

(F.4) 



Using the above, 



isDo 



2/3^2(1 + r7o) /3ao^ 



oo oo 



+ ^ 5Z 5Z ^"^^-"^^ + « 5Z 5Z ^^^^-'-^ 



2/3 



1 



/3 



^/O r 



2/30-2(1 +r7o) /3ao 



^' 



^ V(-)' + - — V - + — V A«e(°) 



^ ^ n=l ^ ^ r>0 ^ ^ ra=l ^ ^ ra=l 






/3cro 



n=l 



/3cro 



r>0 



/3o-o 



E(^^^ 



)^2 



r>0 



- — ^ sign(£) ^ sign(£ - n)sign( 



n] 



£=^=0 



n=/=0,i 



{£ - n)n 27r 



9^ v^ sign(£) -^^ . 

— 2^ — - — 2^ sign(r)sign(£ - 



+ ^ E E A« Ajief + ^ E E Ai°) Afj„e« + 1 E E si'^arXAf 



/3 



2/3 



2 



+ ^EE^^°^^tlA« + 0(/3- 



^7^0 r 



-9/5N 



/3 



£7^0 r 



-« E4:^f-E'ii' -« 



vn^^O 



An 

J 



\n=l r>0 / 



(F.5) 
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We see that the first two terms above cancel with the second and third term in PFq. 

G High Temperature Analysis of One-Loop Truncated 
J\f=2 Gauge-Fixed BFSS 

In this appendix we will present the high temperature expansion to the subleading order of 
the M = 2 Schwinger-Dyson equations in Appendix |B. 4 



G.l Subleading Corrections to the Self-Energies 

For the zero modes, the next-to-leading correction to the self-energies satisfies the following 



equations (see Appendix B.4), with a,b,c,d defined in (6.2), 



(1)M _ 4 V^/_^\2 _ £ V^/ P \2 I ^ V^/ /^ \2 1 _ ^ V^/ P n2 



2/x(i) 4 n«^ 12 f U^^^'' S«'^ 






3 62 52 a\d+l) /3i/2a(t;+i)2' 

j-r{l)V ^ 2 >r^. (3 .2 _ ^ Y^f JL^2 I 1£ V^/ /^ ',2 _ 1^ ^f JL^2 

° (3 ^^27im' I3^^27ir' (3 ^^2Txm' (3 ^^vrr^ 

_2^_2ni^_i4n^ ^^-^^ 

/3/i2 /3 (no^)2 /3 (n^o^ 
_ _8/3 _ V^ _ 2n(^^^ _ i4n^^^^ 

~ 3 62 52 ^2 ' 

^(,)M_ i2n^^)^_ i2/3V2nW^ 
"° (3 (n*0^ «^ 

/i(i) _ 4/3 2n(^)^ 7n«^^ 



2 3 62 ^2 

where we have used 



E 



7r2m2 3' 



Z—/ 'rr2q 



(G.2) 



vr^r2 = ^• 
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The solution to these algebraic equations are, 

n«^' = -2.52/3, n^'^^ = -0.32/3, S^'^"^ = 0.68/3^/2, ^« = 5.53/3. (G.3) 

Now for nonzero modes, 
]t(i)m ^ 4 >^ /3 2 4 >^ /3 2 12 >^ /3 2 1 I ''"^ V^ ^ ^ 



/3 f;^;27rm^ /3 ^^27rr^ /3 ;^;27rm^ 1 + S^_„ /3 ^ 27rr 2n{n - r) 

4 2vm >r^ /3 3 /3 _ 2^ , 8^ , Af A^^nM 
P /3 V 2^^ 27r(n-r) /3/i2 ^ /3/i2 ^ /3/i47rn^ "" 



6u(i) 4n^^^^ n^^^^ 

— Rr^ -I- '^ — _ iQ -'-'-o 
"^ " + 52 ^2 ^2 ' 

TT(i)y ^ 2 >^ /3 2 3 >^ /3 2 14 ^ /3 n2 , ^^ y^ /3 /3 



B ^-^ 271171 B ^-^ 2T:r B ^-^ 27rm B ^-^ 2Txr 2Tx{n — r) 



-y 



/3 /3 2/i(i) , 8^(1) ,8^/3 ^^^y 2 Xlf"" 14 n^^^^ 



/3 ^ 27rr 27r(n - r) /3/i2 ' /3;i2 ' /3;x'27rn' "" /? (H^)^ /^ W)^ 



6^(1) 2n(^)^ Mn^"' 



52 52 a2 



71 



/^ ^ 2 >r^. /3 .2 ^ Y^cJL^2 I 1 V^/ ^ Vl _ £ Y^('A_A2 I 4 \p/ /^ \2 

2 13 ^^2TTm' 2(3 ^^2nr' (3 ^^2nm' (3 ^^2nm' (3 ^^2Txr' 



(3^^2-117X1' (3^^2-nr' (3 ^^2-Km' 2/3 ^^vrr 



2 ni^)^ 7 n«^ 



/3 {uyy /3 (n^)2 

1, 2n(^)^ 7n«^_4/3 2n(^)^ rn^^^^ 



where 



3 52 a2 3 52 a 



j\^ 4 2/9 2(a-66) 



2 



" 3 7r2n2 Va2 + 6 ab^ 

" 12 7r2n2 V 62 



(G.4) 



(G.5) 



We have used 



7r2 



El ^-^ 1 71 



r{n — r) ' -^^ r3(n — r) rfi 
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Therefore, from (6.3) and (G.3), 

i2nj^)^ 4n«^ 6/i(i) 



/3-inw^^ 


4 1/18 
~ 3 7r%2 \^a2 + Q + 




= 5.24183 + "'"^f « 


/3-^n«^ 


7 2 / 8\ 
12 ^ TT^n^ V &V 




- 1 ^650- 1 °-'''''' 




n"' 



a62 J pa^ (5h^ ^ pb' 



/3a2 /362 + /352 



,2 



(G.7) 



G.2 Free Energy in the High Temperature Limit 

To leading order in P, the free energy receives only contribution from the bosonic sector. 
Using the fact that 11^ and 11^ scales as P~^'^ (6.4), the same as their corresponding zero 
modes, we have 



I I 



log/i + 21ogsmh ^ "^ — + 71ogsmh ' '^^ " ^ i ^ ^ 



2°' °V2y °V2 

~ const + 6 log p. 

The subleading contributions to the free energy come from the 0{P^^'^) corrections to the zero 
mode self-energies and also the leading nonzero mode self-energies (relative to the kinetic 
energy), 

. X n(^)^ jytWM 7 ^(l)M . (1) 



j-j{o)y 2n(°)*^ 2i + s^°^*^ 2/i(o) 

+ ^E=f'-E(^w.25:A,. 

/3Fi^) = _ y ( A)2nv' + y ( A)3s^v _ ! y (^)2nA^ + 7 y A^f 

^ 2 ^^2nr ' ^^^27rr^ ' 2 ^^^27r/^ ^ ^ 27ir ' 

1^0 r 1^0 r 



(G.9) 






(G.IO) 
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Note that to this order 



Hence 



"0 



1 I ^M 
i -t- --0 



const 



^(l)M 



^(1)M^(0)M 



(G.ii: 



/3F(^^ + /3F(^^ 



^3/2 



-'"'-0 _|_ ' -'"'-0 _|_ ^ 

/36 2 /3a 



"0 " 



+ 



lu« 



2/33/2(1 + ^)2 2 



/3F3 and /3F4 can be similarly computed to be 



m 



(1) 



+ 



+ 



14 
21 



+ -E 



El 



4 

J^ 14 

27rr^ 



•27rr^ 
/3 ^ 



42 



ms 



M\2 



/3n, 



pf^E' 



A)2_ 

27rr^ 



/3/i^^7r/' 
27r/^ 



+ 



E 
-E 



1^0 

1 

2/3"/i 



A)2 



14 

7 

iSfj, 

42 ^,J^^, 
27rr 
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Therefore, 

^^7^(1) ^ (10554 _ 233862 _ 2220) /33/2 

^864^ '' (G.15) 

~ - 3.89872/3^/2. 

Hence 

PP = const + 6 log/3 - 3.89872/3^/2 ^ (^^^3^)_ ^q^^q^ 
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